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ABSTRACT 

The  purpose  of  this  study  was  the  development  of  a  computer  system 
for  research  and  instruction  in  signal  processing  for  underwater  sound 
applications. 

The  nature  of  the  problem  demanded  the  use  of  three  computers  inter- 
faced by  hardware  and  software  to  operate  as  a  single  system.  A  general 
purpose  analog  computer  was  employed  to  generate  signals  and  noise, 
provide  analog  processing  and  filtering  and  to  allow  the  input  of  actual 
signals  collected  in  the  field.  A  general  purpose  digital  computer 
provided  for  the  examination  of  existing  processing  methods  and  for  the 
investigation  of  new  techniques.  A  high  level  graphics  computer  was 
used  for  data  display  and  as  the  system  control  station  which  provided 
for  a  high  degree  of  human  operator  processing  and  interaction. 

All  the  hardware  required  was  provided  in  the  Computer  Laboratory 
of  the  Naval  Postgraduate  School,  Department  of  Electrical  Engineering 
and  this  study  prepared  the  software  for  its  convenient  and  efficient 
use  for  this  application. 
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I.   INTRODUCTION 

A.  PURPOSE 

The  purpose  of  this  work  was  to  develop  a  convenient  method  and 
efficient  programs  to  employ  the  unique  computer  system  of  the  Naval 
Postgraduate  School  Computer  Laboratory  to  build  a  hybrid  system  for 
signal  analysis  and  processing.  The  primary  good  of  the  system  was  to 
achieve  a  high  speed  computational  tool  with  a  high  degree  of  human 
operator  interaction.  In  addition  the  system  was  to  provide  the  flexi- 
bility and  interaction  which  would  permit  the  operator  to- use  his 
heuristic  capabilities  in  approximating  solutions  in  a  design  procedure 
when  the  search  for  adequate  patterns  and  processors  is  the  goal,  or  to 
permit  the  student  to  verify  theoretical  results  for  a  wide  variety  of 
situations  that  he  may  configure  almost  at  will. 

It  was  also  assumed  from  the  beginning  that  the  system  should  be 
directly  aimed  to  the  underwater  sound  detection  problem.  This  decision 
not  only  dictated  the  name  of  the  system  namely  a  Hybrid  Sonar  Simulation, 
but  also  directed  the  choice  of  the  preprogramed  applications  initially 
included.  This  choice  does  not,  however,  restrict  the  application  of  the 
program  to  a  much  wider  field  of  signal  processing  limited  primarily  only 
by  the  user's  imagination  and  expertise. 

B.  SYSTEM  REQUIREMENTS 

After  establishing  the  general  purpose,  some  specifications  were  set 
up  to  further  specify  the  scope  to  be  achieved. 

In  order  to  enable' the  user  an  easy  and  rapid  perception  of  the 
signals  and  the  results  of  their  analysis  or  processing,  a  graphical 


display  was  chosen.  To  allow  him  to  interact  with  the  process,  it  was 
necessary  to  make  available  means  of  controlling  the  signals,  processors 
and  operations  at  the  same  location  where  the  situation  was  being  dis- 
played and  while  it  was  shown.  The  fundamental  requirement  for  this  type 
of  operation  is  very  high  speed  computation. 

Since  comparison  of  signals  and  patterns  was  one  of  the  main  objec- 
tives, the  user  should  be  able  to  choose  among  different  classes  of 
signals  and  then  vary  their  parameters.  Also  the  characteristics  of 
filters  should  be  under  control. 

For  analysis  and  processing  there  was  established  the  requirement 
for  spectral  analysis,  auto-  and  cross-correlation  and  convolution. 

Three  types  of  processors,  a  linear  cross-correlator  and  two  non- 
linear processors  were  included  as  well  as  the  capabilities  to  determine 
signal-to-noise  ratio,  input-output  performance  characteristics  ("proces- 
sing gain")  and  the  receiver  operating  characteristics  (ROC). 

Random  noise  was  necessary  to  be  mixed  with  the  uncorrupted  signals 
in  controlled  amounts.  A  Gaussian-like  noise  was  considered  adequate 
for  initial  test  purposes. 

Finally  a  possibility  for  taking  hard  copies  of  the  curves  displayed 
had  to  be  incorporated  to  allow  record  keeping  for  further  studies  and 
comparisons. 

C.  ELEMENTS  AVAILABLE 

The  elements  available  in  the  computer  laboratory  and  used  to  compose 
the  system  include  (see  Fig.  1.1): 

Digital  computer  -  XDS  9300  -  a  medium  sized  general  purpose  computer 
with  32K  words  of  main  memory  and  a  processor  capable  of  handling  a  large 
set  of  operations  with  very  good  efficiency.  In  the  installation  it  is 
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connected  to  the  analog  computer  and  the  graphics  digital  computer  which 
allows  the  hybrid  configuration.  Also  available  are  several  of  the 
standard  input-output  devices  among  them  a  magnetic  drum  secondary 
storage  with  a  two  million  character  capacity. 

Analog  computer  -The  Comcor  CI  5000  is  a  very  well  equipped  analog 
computer  which  includes  a  parallel  logic  portion  that  allows  extensive 
decision  making  and  control  functions  to  be  performed. 

Graphics  computer  —  The  AGT-10  is  a  small  general  purpose  computer, 
with  8K  words  main  memory  which  has  as  a  special  feature  a  CRT  output  and 
the  necessary  hardware  to  write  on  the  screen,  generate  vectors  and 
refresh  the  pictures. 

Converters  —  Both  a  high  precision  analog  to  digital  converter,  or 
sampler,  and  digital  to  analog  converters  are  available. 

Software  —The  system  software  includes  the  normal  mathematical  func- 
tions and  subroutines  to  perform  some  FORTRAN  operations,  a  set  of  programs 
to  control  the  analog  computer  from  the  digital,  programs  to  control  the 
graphics  terminal,  two  subroutines  to  plot  curves  with  the  line  printer 
and  a  program  to  compute  fast  Fourier  transforms  according  to  the  so  called 
Cooley-Tuckey  method  [Ref.  2]. 

D.  GENERAL  APPROACH  TO  THE  SOLUTION 

The  logical  approach  was  to  use  the  elements  available  according  to 
their  best  performance  characteristics  as  applied  to  the  problem. 

Naturally,  the  graphics  computer  was  taken  as  the  data  presentation 
component  and  then  became  also  the  system  control  unit.  Sitting  in  front 
of  it  the  engineer  incorporates  himself  into  the  system  as  the  center  of 
decisions.  He  selects  an  operation,  analyzes  the  results,  determines 


the  next  step,  compares  the  new  results,  and  introduces  new  parameters 
until  satisfied  with  the  information  accumulated.  Since  the  graphics 
terminal  is  itself  a  general  purpose  computer,  the  types  of  data  to  be 
displayed  as  well  as  the  format  of  the  display  may  be  varied  to  investi- 
gate the  effect  on  Human  Operator  performance. 

In  nature  the  great  majority  of  data  or  information  comes  in  a 
continuous  or  analog  form,  the  discretization  being  an  human  abstraction 
or  simplification  for  easier  handling.  In  particular,  that  is  the  situ- 
ation encountered  in  the  detection  problem  where  the  signals  are 
continuous.  For  these  reasons  the  analog  computer  was  utilized  to 
generate  the  signals.  The  analog  computer  offers  the  possibility  of 
wider  range  of  variations  with  very  simple  and  easy  to  perform  change  of 
parameters.  By  using  signals  of  the  same  nature  as  the  real  ones, 
practically  no  changes  are  required  if  instead  of  simulated  functions 
true  signals  are  the  object  of  analysis  either  in  the  form  of  a  pre- 
recorded reproduction  or  as  the  output  of  a  receiver. 

To  analyze  or  process  the  signals  the  discretization  becomes  a 
necessity  since  the  digital  machines  are  by  far  more  precise  and  wery 
well  adapted  to  the  computations  described  by  fixed  algorithms.  Hence 
the  choice  of  the  digital  computer  for  this  part  of  the  system. 

Two  ways  may  be  used  to  show  the  arrangement.  From  a  functional  stand- 
point the  operations  and  their  inter-relationship  are  the  constituent 
blocks  as  presented  in  Fig,  1.2.  Depending  on  the  manner  these  functions 
are  employed,  two  modes  of  operation  result  and  may  be  named  ANALYSIS  and 
SIMULATION.  In  the  ANALYSIS  mode  it  was  intended  that  basically  one  func- 
tion be  performed  at  a  time  under  operator  control  thereby  giving  him  an 
opportunity  to  verify  intermediate  results  and  to  proceed  accordingly. 
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Through  the  SIMULATION  mode  some  pre-programmed  problems  were  configured 
and  automatic  sequences  could  be  entered  to  determine  processing  gain, 
receiver  operating  characteristics  and  detection  distance.  Some  of  the 
processing  functions  could  also  be  executed  by  a  program  that  performs 
correlation  by  a  point  by  point  multiplication  and  addition;  since  these 
operations  serve  no  other  purpose  than  to  demonstrate  the  long  time  re- 
quired by  this  direct  approach  they  were  thought  of  as  belonging  to  a 
DEMONSTRATION  mode. 

From  a  physical  point  of  view  the  system  is  described  as  depicted  in 
Fig.  1.3.  The  digital  computer,  besides  providing  most  of  the  processing 
functions,  exerted  effective  control  of  the  other  units  from  which  it 
received  interrupts  to  initiate  or  synchronize  operations. 

E.  ELEMENTS  CREATED 

To  implement  this  system  several  elements  were  created  to  execute  the 
functions  and  satisfy  the  requirements  established.  A  brief  description 
of  these  elements  is  presented  in  the  following  paragraphs. 

1 .  Control 

A  program  was  written  to  read  the  control  codes  typed  on  the 
graphics  terminal  and  compare  them  against  a  table  to  establish  pointers 
that  indicated  the  proper  sequence  to  be  executed. 

Another  program  was  introduced  to  allow  the  operator  full  control 
of  the  analog  computer  from  the  graphics  terminal.  Through  this  feature 
it  is  possible  to  select  modes,  set  or  test  lines,  set  potentiometers  or 
scan  the  various  analog  components  by  typing  simple  messages  at  the 
console, 

2.  Signal  Generation 

Signal  generation  was  developed  from  sine  generators  patched  on 
the  analog  board.  Bearing  in  mind  the  sonar  application,  one  of  these 
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generators  was  associated  with  additional  analog  amd  logic  circuitry  to 
permit  simulation  of  CW,  FM-slide  or  FSK  codes.  This  generator  was  also 
equipped  with  positive  feedback  to  make  possible  a  more  realistic  effec- 
tive sampling  rate  which  required  the  analog  computer  to  run  at  a  very 
high  frequency.  Without  some  regeneration  in  this  case  the  wave  damping 
would  be  intolerable. 

Three  other  sine  generators  were  implemented  to  simulate  propeller 
signatures,  which  present  a  low-pass  line  spectra  composed  of  a  funda- 
mental and  several  of  its  harmonics  [Ref.  1].  To  achieve  this  result  the 
sine  waves  were  fed  to  a  combination  of  multipliers  whose  outputs  were 
summed  in  controlled  amounts  to  give  a  final  pattern  with  up  to  six 
harmonics,  each  of  which  with  any  desired  intensity. 

3.  Noise  Generation 

To  simulate  a  random-like  noise  a  pseudo-random  binary  sequence 
was  generated  by  a  maximal  length  feedback-shift  register  and  then  trans- 
formed into  a  pseudo-random  analog  square  wave.  This  signal,  after  proper 
filtering  presents  Gaussian  statistics  at  least  up  to  a  first-order  [Ref.  6] 

4.  Acquisition 

Signal   acquisition  to  the  processor  was  normally  made  after  filter- 
ing and  then  through  the  sampler  or  Analog-to-Digital   converter  (ADK). 

Two  filters  were  included:     an  RC  lowpass  aand  a  narrowband,  high-Q, 
second-order  filter.     The  latter  offered  the  advantage  of  being  easily 
controlled  from  the  graphics  terminal   to  adjust  frequency,  bandwidth  or 
gain. 

The  purpose  of  the  filters  was  to  represent  transducer  and  part 
of  preprocessor  characteristics,  but  they  could  also  simulate  or  be  pre- 
ceded by  other  filters  with  any  other  objective,  such  as  transmission 
medium  characteristics. 
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The  sampler  was  controlled  by  a  program  that  not  only  synchronized 
the  operation  but  also  simulated  the  remaining  parts  of  a  pre-processor 
for  the  non-linear  cases,  and  stored  the  samples,  conveniently  normalized, 
into  the  proper  locations  of  the  processor  memory. 

The  samples  were  distinguished  and  stored  as  belonging  either 
to  the  REPLICA  or  to  the  SIGNAL.  While  the  choice  of  names  relects  the 
application  to  sonar,  for  a  general  case  these  names  serve  only  to  identify 
independent  storage  areas  where  signals  could  be  stored  for  further  opera- 
tion or  retrieval.  Similarly,  a  third  area  where  the  results  of  correlations 
were  placed,  received  the  name  of  CORR. 

5.  Spectral  Analysis 

A  special   subroutine  was  prepared  to  make  use  of  the  Fast  Fourier 
transform  program  and  prepare  the  results  for  display. 

6.  Processing 

A  replica  correlator  simulation  was  employed.  Depending  on  how 
the  signals  have  been  preprocessed  the  processor  became  a  linear, 
amplitude-limited  or  a  clipper-limited  type. 

The  program  written  to  execute  the  correlation  was  based  on  the 
properties  of  discrete  Fourier  transforms  and  on  the  time  domain  convolu- 
tion theorem  as  applied  to  the  discrete  case.  For  cross-correlating  the 
replica  against  a  very  long  time  series,  as  is  the  case  where  processing 
the  incoming  signal  in  a  detection  problem,  a  partition  technique  was 
adopted  [Ref.  4]. 

The  correlation  could  also  be  evaluated  by  multiplying  both 
functions,  adding  the  results,  shifting  and  repeating  these  operations 
until  the  replica  had  been  slid  entirely  across  the  signal.  The  time 
required  for  this  procedure  was  incompatible  with  the  real-time  specifi- 
cation and  the  method  was  left  for  illustration  only. 
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7.  Performance  Determination 

Under  this  general  classification  were  considered  the  input-output 
signal-to-noise  ratio  characteristic,  the  receiver  operating  characteris- 
tic and  distance  measurement. 

A  program  was  prepared  to  automatically  sample  the  replica  and 
the  signal,  form  the  signal  either  with  noise  alone  or  with  noise  plus 
transmission  modulation  code,  process  the  received  samples,  compute  power, 
determine  threshold  and  threshold  crossing,  evaluate  probability  of 
detection  or  false-alarm  and  measure  distance,  all  in  the  proper  sequence 
according  to  the  case. 

Power  determination,  when  necessary,  was  computed  through  the 
application  of  the  Wiener-Kinchin  relation  and  Parseval rs  theorem. 

The  threshold  was  normalized  to  the  peak  value  of  the  correlator 
output  when  no  noise  was  present. 

When  using  these  automatic  operations,  where  only  the  final 
results  and  not  the  waveforms  were  of  interest,  the  sampling  rate  utilized 
was  related  to  twice  the  information  bandwidth.  This  allowed  processing 
samples  with  a  more  realistic  information  content.  Attempts  to  display 
the  waveforms  at  these  sampling  rates  showed  that  while  the  signals  were 
completely  distorted,  they  revealed  amplitude  variations  at  a  rate 
comparable  to  the  significant  bandpass.  For  signal  analysis,  when  the 
waveforms  were  required  to  be  displayed  without  great  distortion,  the 
sampling  rate  had  to  be  considerably  greater  than  twice  the  highest 
frequency  present. 

8.  Display 

Programs  were  prepared  to  display  the  functions  and  results 
covering  all  cases  treated. 
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Figure  1.4.  Spectrum  of  a  CW  signal.  This  figure  is  a  photograph  of 
the  actual  display.  It  shows  the  time  function,  absolute  value  and 
phase  angle  of  the  transform.  The  location  of  the  impulses,  at  15000  Hz, 
serve  as  reference  for  the  frequency  plots  in  the  next  pictures. 
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Figure  1.5.  Spectrum  of  an  FM  signal  formed  by  a  5000  Hz  sinu- 
soidal linearly  modulated  by  1000  Hz.  Comparison  with  Fig.  1.4 
shows  the  wider  band  covered.  Discrepancies  observed  in  the  phase 
angle  plot  here,  as  in  some  of  the  following  pictures,  are  due  to 
digital  computer  noise  in  the  computation  of  arctan. 
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Figure  1,6.  Spectrum  of  a  sinusoidal  signal  of  5000  Hz  two^tone 
frequency  modulated  by  t500  Hz. 
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Figure  1.7.     Spectrum  of  unfiltered  noise.     The  unfiltered  noise 
is  a  pseudo-random  square  wave  which  has  been  passed  only  through 
a  low-pass  filter  with  cut-off  frequency  beyond  the  region  of 
interest.     This  filtering  is  necessary  to  assure  Gaussian  first- 
order  statistics  to  the  simulated  noise.     In  the  band  up  to  5000  Hz 
Ccompare  with  Fig.   1.4)  no  significant  peaks  appear  in  the  spectrum. 
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9.  Copying 

A  subroutine  was  written  to  copy  the  data  displayed  on  the  screen 
into  the  format  required  to  use  the  two  programs  for  line  printer  plotting. 

Also  a  sequence  in  this  subroutine  was  prepared  to  send  the 
information  to  the  analog  computer,  through  the  digital  to  analog 
converter,  in  order  to  plot  the  curves  with  a  paper  recorder.  This  tech- 
nique proved  very   useful  in  several  applications, 

F.  SYSTEM  CAPABILITIES 

The  system  was  made  flexible  enough  to  allow  a  wide  range  of  appli- 
cations. Some  of  these  were  checked  with  very   impressive  results  as 
demonstrated  below. 

Initially,  Figures  1.4  through  1.8  show  types  of  signals  used. 
Fig.  1.4  presents  a  CW  burst  of  5000  Hz  and  2.53  msec;  also  shown  is  its 
spectrum,  magnitude  and  phase.  In  Fig.  1.5  the  5  kHz  signal  is  linearly 
FM  modulated  by  1000  Hz  and  in  Fig.  1.6  the  5  kHz  carrier  has  an  FSK 
modulation  of  ±500  Hz,  Fig.  1.7  depicts  a  2.53  msec  sample  of  unfiltered 
noise;  the  frequency  plot  covers  the  range  of  ±50  kHz,  Fig.  1.8  is  a 
simulated  propeller  signature  with  a  fundamental  of  60  Hz  and  six  harmonics 

Figure  1.9  shows  an  example  of  filter  characteristics  determination 
through  the  application  of  an  impulse  at  its  input.  In  the  upper  part  the 
impulse  response  is  presented  followed  by  the  magnitude  and  phase  of  the 
transfer  function.  The  filter  tested  was  a  narrowband  centered  at  5  kHz 
with  a  bandwidth  of  less  than  1  kHz. 

In  the  next  figure,  Fig.  1.10,  the  same  random  noise  previously  dis- 
played in  Fig.  1.7  is  seen  but  after  being  passed  through  the  narrow 
bandpass  filter.  This  experiment  was  a  demonstration  that  noise  behaves 
as  a  sinusoidal  signal  when  narrow-band  filtered. 
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Figure  1.8,  Representation  of  a  low^pass  periodic  signal  similar 

in  characteristics  to  a  propel ler^generated  noise.  The  line  spectrum 

shows  well  defined  harmonics  up  to  the  sixth.  The  fundamental 
frequency  was  60  Hz. 
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Figure  1.9.     This  figure  presents  the  narrow-band  filter  impulse 
response,  and  both  the  amplitude  and  phase  angle  of  its  transform 
or  transfer  function. 


23 


Figure  1.10.  In  this  picture  the  noise  spectrum  is  seen  again  but 
now  after  narrow-band  filtering.  The  spectral  density  is  similar 
to  that  of  a  sinusoidal  signal  but  the  phase  angle  has  a  random 
variation. 
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Figure  1.11.     The  FM-signal  of  Fig.   1.5  is  shown  after  passed 
through  the  narrow-^band  filter.     The  effect  of  the  filter 
narrow-pass  is  reflected  in  the  time  function  by  a  reduction 
in  amplitude  for  the  frequencies  above  and  below  the  center 
frequency.     In  the  frequency  domain  a  reduction  in  bandwidth  is 
observed.     The  sample  in  this  example  is  twice  as  large  as  that 
of  Fig.   1.5  and  shows  the  effect  of  the  transient  filter  response, 
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Figure  1.12.     This  picture  shows  how  the  filter  output  of  the 
previous  figure  can  be  obtained  by  computation.     Above  is  displayed 
the  filter  impulse  response  time-reversed.     In  the  middle  the 
unfiltered  FM-signal,  and  below  the  result  of  the  convolution  of 
these  two  functions. 
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Figure  1,13.     Power  spectral  density  of  the  signal   presented  in 
Fig.   1.8  but  now  highly  contaminated  by  noise,     The  auto-correlation 
function  and  its  transform  are  shown. 
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Figure  1.14.     Auto-correlation  function  of  a  CW-pulse. 
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Figure  1.15.     Auto-correlation  of  the  FM  signal   shown  previously 
in  Fig.   1.5.     Comparison  with  the  auto-correlation  of  the  CW 
signal   presented  in  Fig,   1.15  shows  the  better  definition  of  the 
correlation  as  a  consequence  of  the  wider  bandwidth. 
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Figure  1.16.  Autocorrelation  of  the  FSK  pulse  presented  tn  Fig,  1.6, 
By  comparison  with  Fig.  1.15  it  is  seen  that  the  uneven  distribution 
of  the  energy  over  the  bandwidth  affects  the  result  of  the  correlation, 
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Figure  1.17.     Autocorrelation  of  unfiltered  noise.     The  picture 
clearly  shows  the  lack  of  correlation  among  random-noise  samples. 
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Figure  1.18.     In  this  picture,  the  CW  signal  of  Fig.   1,4  was 
amplified  and  clipped  after  having  been  passed  through  a  narrow- 
band filter.     The  envelope  of  the  auto-correlation  is  similar  to 
that  of  the  undistorted  signal,  which  shows  that  narrow-band 
signals  may  be  treated  by  absolute-value  processor  models. 
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In  Fig.  1.11  the  FM  signal  of  Fig.   1.5  is  again  displayed  but  now 
after  filtering.     In  Fig,  1.12  the  same  result  is  obtained  by  convolving 
the  filter  impulse  response  with  the  unfiltered  signal  —  the  upper  part 
of  the  figure  is  the  filter  response  flipped,  the  central   the  unfiltered 
signal  and  the  lower  the  result  of  convolution.     This  experiment  demon- 
strated the  capability  of  the  system  to  simulate  digital   filters  or  any 
other  filter  once  the  impulse  response  is  known. 

Figure  1.13  shows  the  result  of  power  spectral  density  evaluation  of 
the  propeller  signature  already  mentioned  but  now  noise  added  and  after 
lowpass  filtering.     Displayed  are  the  auto-correlation  of  the  signal   and 
the  magnitude  and  phase  of  the  power  spectrum. 

The  next  four  figures,  1.14,  1.15,  1,16  and  1.17  depict  the  auto- 
correlation of  the  signals  previously  shown  in  Figures  1.4  through  1.7, 
i.e.,  CW,  FM,  FSK  and  noise. 

Figure  1.18  shows  the  auto-correlation  of  a  CW  signal   processed  by  a 
clipper-limited  type  of  processor.     Since  only  the  zero  crossing  informa- 
tion is  kept,  the  resultant  signal   is  a  binary  signal. 

The  following  figures  present  the  results  of  the  cross-correlation 
of  a  replica  and  a  signal   composed  of  the  transmission  pattern  plus  noise, 
The  cross-correlation  was  normalized  in  such  a  way  that  if  the  transmis- 
sion code  only  was  present  the  peak  output  would  be  one.     All   samples  are 
taken  after  narrowband  filtering  and  correspond  to  the  FM-slide  signal. 

Figures  1.19  and  1.20  refer  to  a  linear  cross-correlator  output  for 
two  different  signal-to-noise  ratios  while  Figs.   1.21   and  1.22  use  the 
same  SNR  but  applied  to  a  clipper-limited  processor.     The  experiment 
showed  that  for  high  signal-to-noise  ratio  the  non-linear  processor  has 
a  much  lower  sensitivity  than  the  linear  one;  the  difference,  however, 
is  almost  non-existent  for  a  low-input  signal-to-noise  ratio. 
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Figure  1.19.     Cross-correlation  of  an  FM-replica  against  this 
signal   plus  noise  at  a  high  signal-to-noise  ratio.     The  periodic 
or  circular  nature  of  the  Discrete  Fourier  Transform  is  evident 
from  the  result  of  the  correlation. 
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FiQure  1  20.  This  picture  is  the  result  of  the  same  operation  as 
that  used  to  obtain  Fig.  1.19.  but  applied  when  the  signal -to-noise 
ratio  is  low.  Although  the  functions  are  normalized  to  the  same 
maximum  amplitude  in  the  display,  the  deterioration  caused  by 
noise  is  evident,  especially  when  the  central  "bulge"  is  observed. 
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Figure  1.21,  This  figure  shows  the  result  of  processing  the 
signals  of  Fig.  1.19  with  a  clipper^limi ted  processor. 


36 


Figure  1,22.     This  figure  presents  the  results  of  applying  a 
clipper-limited  processor  to  the  signals  of  Fig.   1.20,  where  a 
low  signal-to-noise  ratio  was  employed.     Since  the  amplitudes 
are  displayed  normalized  to  the  same  value,  comparison  can  be 
made  by  considering  the  peak-to-peak  oscillations  of  the  outer 
segments.     At  low  SNR  the  non-linear  processor  approaches  the 
linear  processor  better  than  when  the  SNR  is  high. 


37 


For  the  next  experiment  CFig.  1.23)  in  order  to  simulate  a  doppler- 
shifted  reception,  the  FM  signal  was  modified,  but  not  the  replica.  The 
input  signal-to-noise  ratio  is  the  same  as  for  the  case  of  Fig.  1.19  in 
order  to  present  the  degradation  caused  by  doppler. 

In  Fig.  1.24  a  replica  of  a  propeller  signature  is  cross-correlated 
against  a  long  sequence  of  incoming  samples  where  the  desired  signal  is 
present  but  at  a  very  low  signal-to-noise  ratio.  The  replica  was  corre- 
lated against  the  long  sequence  using  the  partition  method;  the  result 
displayed  is  only  a  section  of  the  correlation  corresponding  to  9.12 
msec.  It  was  observed  that  the  peak  output  was  low,  although  the  oscil- 
latory character  of  the  correlation  seemed  to  carry  much  of  the  information 
required  for  detection.  A  Fourier  transform  of  the  result,  illustrated 
in  Fig.  1.25,  showed  how  this  information  could  be  identified. 

The  last  figure  of  this  series,  Fig.  1.26  is  the  result  of  an  experi- 
ment similar  to  that  of  Fig.  1.20  but  with  the  sampling  rate  just  greater 
than  twice  the  bandpass  (up  to  this  point  the  sampling  rate  used  was 
twenty  times  the  highest  frequency).  It  was  noted  that  the  displayed 
replica  is  completely  distorted  (no  narrow  band  filters  were  used  to 
recover  the  signal  from  the  samples)  but  revealing  an  amplitude  variation 
in  the  order  of  the  bandwidth.  The  cross-correlation,  although  of  unrecog- 
nizable shape,  has  a  highest  peak.  In  this  case,  where  both  the  total 
number  of  samples  and  the  bandwidth  have  been  conserved,  the  reduction  in 
the  sampling  rate  corresponded  to  a  larger  integration  time  with  less 
redundant  samples  being  processed.  Conversely,  in  the  case  of  Fig.  1.20 
the  integration  time  was  short  and  several  non-independent  samples  were 
processed,  contributing  nothing  to  the  result. 

The  importance  of  input-output  performance  and  ROC  determination  for 
comparison  of  processors  and  transmission  codes  was  verified  by  several 

problems  simulated  as  described  in  a  later  section, 
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Figure  1.23.  This  figure  presents  the  result  of  an  experiment 
designed  to  simulate  a  doppler  shift  of  the  incoming  signal. 
Comparison  with  Fig.  1.19  demonstrates  the  degradation  in  the 
correlation  as  a  consequence  of  the  shift. 
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Figure  1.24.  In  this  picture  a  section  of  the  cross-correlation 
of  the  propeller  noise  replica  and  the  continuously  received 
signal  with  large  amount  of  noise  is  shown.  Note  that  the 
cross-correlation  function,  if  not  further  processed,  does  not 
reveal  any  information. 
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Figure  1.25.     This  figure  shows  the  result  of  the  spectral  analysis 
of  the  cross-correlation  function  displayed  in  Fig.   1.24. 
Comparison  with  Fig.  1.8  demonstrates  that  the  presence  of  the 
original   signal   is  now  readily  identifiable.     The  same  is  not  true, 
however,  in  the  power  spectrum  of  the  signal  before  correlation 
(Fig.   1.13). 
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Figure  1.26,  This  figure  presents  the  effect  of  sampling  the 
signals  of  Fig.  1.20  at  a  sampling  rate  close  to  twice  the  signal 
bandwidth.  Note  that  while  the  signal  is  badly  distorted  at  this 
low  sampling  rate,  and  without  any  reconstruction  filtering,  the 
peak  of  the  cross-correlation  is  much  better  defined. 
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G.  DETAILED  DESCRIPTION 

In  the  following  sections  a  detailed  description  of  the  system  imple- 
mentation together  with  a  review  of  the  mathematical  concepts  utilized 
will  be  given.  Further,  a  brief  section  on  the  graphics  terminal  and 
discussion  of  sample  simulations  is  presented.  Finally,  some  conclusions 
and  recommendations  for  additional  applications  are  provided. 

The  author  does  not  claim  originality  for  the  theory  discussed,  which 
can  be  found  in  good  reference  books  such  as  those  mentioned  in  the  text. 
Perhaps  only  a  certain  degree  of  innovation  exists  in  the  interpretation 
given  to  information  bandwidth  which  seems  to  depart  somehow  from  the 
conventional . 

The  work  described  here  was,  and  was  intended  to  be,  an  engineering 
study.  A  program  requiring  a  solid  theoretical  background,  a  clear  view 
of  the  objectives,  and  a  willingness  to  make  the  unavoidable  compromises 
characteristic  of  real  problems,  in  order  to  transform  the  resources 
available  into  a  useful  product. 
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II.  DIGITAL  PROCESSING 

The  system  was  implemented  to  process  the  signals  digitally,  since 
the  great  majority  of  processors  of  interest  are  digital.  The  reason  for 
that  rests  mainly  upon  the  ability  of  digital  machines  to  handle  algorith- 
mic operations  wery   rapidly  and  with  high  precision.  The  discretization 
of  data  required  for  digital  processing  does  not  necessarily  represent  a 
reduction  in  information.  If  the  sample  rate  is  properly  adjusted,  only 
independent  samples  may  be  processed  and  no  time  is  wasted  with  the 
redundancy  existing  in  the  original  signal. 

To  use  the  XDS-9300  of  the  laboratory  it  was  necessary  to  write 
programs  in  order  to  simulate  the  system  and  processors.  Since  this  is 
a  general  purpose  computer  several  other  processors  can  be  simulated 
without  any  hardware  change,  requiring  only  the  addition  of  new  programs 
or  modification  of  those  initially  included. 

In  this  section  a  description  of  the  computer  program,  its  usage,  and 
the  theoretical  reasons  determining  its  construction  are  presented  together, 
in  order  to  make  it  easier  for  the  reader  to  follow  the  discussion, 

A.  OVERLAY  STRUCTURE 

Since,  as  noted  before,  the  computer  has  to  store  in  core  a  large 

collection  of  library  subroutines  to  be  used  in  the  program,  including 

those  to  control  the  analog  computer  and  the  graphics  terminal,  and  also 

has  to  store  relatively  large  arrays  for  the  computations  of  correlations 
and  convolutions,  little  space  is  left  to  the  program  itself.  To  overcome 

this  difficulty  the  program  was  segmented  as  shown  in  Fig.  2.1. 


44 


UJ 
cr 

O 

a: 
c/> 

>- 
< 

-J 

UJ 

> 
o 


CM 
UJ 

cc 
o 


45 


In  brief,  the  purpose  of  each  segment  is  as  follows.  Segment  MAIN 
contains  the  program  initialization  and  the  pivot  points  where  the  status 
is  saved  when  a  new  overlay  is  called  in.  Segment  G  is  composed  of  the 
subroutines  required  to  control  the  analog  computer  from  the  graphics 
terminal.  Segment  A  permits  bringing  in  functions  and  performing  opera- 
tions on  them  under  operator  control,  and  constitutes  the  base  of  the 
ANALYSIS  mode.  Segment  AB  is  the  principal  part  of  the  processor  and 
computes  the  cross-correlation  of  two  functions;  it  may  be  called  either 
from  ANALYSIS  or  PROCESS  through  a  bridge  in  the  main  segment.  Segment 
P  allows  the  representation  of  a  processor  to  be  automatically  put  to 
work  and  enters  what  is  called  the  SIMULATION  mode.  Segment  T  is  designed 
to  copy  what  is  presented  on  the  screen.  Segment  D  computes  correlations 
through  a  direct  approach  for  a  DEMONSTRATION  mode.  Finally,  segment  MA 
provides  graph  and  test  initialization,  decodes  the  instructions  typed 
in  by  the  operator,  and  selects  the  path  to  be  followed.  This  segment 
provides  the  capability  to  select  the  type  of  processor  or  pre-defined 
signals,  and  to  enter  new  parameter  values. 

If  the  system  is  to  be  modified  by  addition  of  new  segments,  the 
maximum  size  allowed  is  approximately  4500  (octal)  words  per  segment. 
This  limit  however  depends  on  the  addition  or  deletion  of  library  sub- 
routines and  also  on  the  size  of  the  sample  as  will  be  noted  later. 

B.  MAIN  PROGRAM 

The  function  of  the  main  program,  besides  serving  as  a  bridge  between 
the  various  operations  and  providing  initialization  of  the  simulation,  is 
to  define  permanent  location  in  memory  for  some  constants  and  the  arrays 
where  the  samples  and  results  of  operations  are  stored. 
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The  constants  are  stored  in  the  COMMON  area  and  their  meaning  will 
be  explained  with  the  subroutines  in  the  Appendix.  One  exception  to  this 
rule  is  the  constant  M  that  defines  the  sample  size;  it  can  be  modified 
by  the  operator  and  is  initially  set  to  256,  its  maximum  allowable  value. 
The  value  of  M  must  be  a  power  of  2  due  to  the  methods  used  for  computation 
of  correlation  and  Fourier  transforms. 

The  arrays  defined  in  the  main  program  are  passed  to  the  subroutines 
as  arguments;  with  this  approach  it  is  readily  possible  to  make  use  of 
the  same  storage  area  for  different  applications  according  to  the  situ- 
ation, and  it  is  not  necessary  to  modify  the  subroutines  when  a  new  value 
is  used. 

The  sample  size,  M,  serves  as  a  parameter  to  define  the  blocks  of 
array  storage.  At  this  point  it  is  convenient  to  remember  that  in  the 
XDS-9300  a  real  variable  uses  two  words  of  storage  while  an  integer 
requires  only  one.  It  follows  then,  that  one  block  of  M  locations  actually 
requires  2M  words. 

The  map  of  Fig.  2.2  shows  the  correspondence  of  designations  through- 
out the  program. 

The  basic  assignment  of  the  arrays  is  as  follows: 

YR  —  replica,  or  the  signal  transmitted;  size  M 

YS  —signal,  or  received  signal;  size  3M,  what  allows  a  faster 
computation  of  recursive  correlation 

YC  —  correlation  or  result;  size  2M 

VAGO  -  empty  buffer  used  for  several  applications;  size  5M 

IB  —  buffer  initially  designed  for  text  and  graph  blocks;  size  2M 
(plus  one  word) 

Before  any  sample  is  entered,  the  main  program  creates  a  dummy  replica 

and  signal,  two  different  size  rectangular  pulses  which  allow  easy 

verification  of  the  system  functioning.  These  two  functions  are  immediately 
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dumped  into  secondary  storage.  This  is  also  always  done  with  the  replica, 
signal  and  correlation  in  the  ANALYSIS  mode.  However,  only  the  yery   last 
one  of  each  of  these  functions  entered  into  the  system  is  saved. 

C.  SUBROUTINES  SEGSAMP,  XEGSAMP  AND  SAMPLE 

Subroutines  SEGSAMP  and  XEGSAMP  have  the  responsibility  of  controlling 
the  sampling  of  the  analog  computer.  They  are  identical  but  have  distinct 
identifiers  since  they  appear  in  two  different  levels  of  segmentation. 
They  also  provide  normalization  so  that  the  amplitude  of  the  samples  is 
less  than  a  maximum  value  compatible  with  the  system. 

Another  function  of  SEGSAMP  is  to  prepare  the  signals  for  the  simula- 
tion of  three  types  of  processors:  linear  cross-correlator,  amplitude- 
limited  cross-correlator  and  clipper-limited  cross-correlator. 

SEGSAMP  controls  three  types  of  sampling  according  to  the  value  of 
INDEX,  i.e.,  replica,  initial  signal,  or  continuation  signal. 

The  replica  normalization  factor,  RCOF,  signal  normalization  factor, 
SCOF,  and  their  ratio  Q  are  saved  in  common  storage,  because  they  are 
necessary  to  the  computation  of  the  input  signal-to-noise  ratio,  SNR; 
Q  is  also  necessary  to  correct  the  computed  threshold  value  for  detection. 

The  processor  type  is  determined  according  to  ISAMP,  i.e.,  linear, 
clipper-limited  or  amplitude-limited. 

INDEX  has  the  initial  value  of  one  and  ISAMP  is  zero;  these  numbers 
however  may  be  set  as  desired  by  the  operator. 

The  actual  samples  are  not  taken  by  SEGSAMP.  After  the  subroutine 
is  entered  and  the  proper  sampling  sequence  is  selected,  the  computer 
is  trapped  in  a  loop  until  an  interrupt  is  generated  by  the  sample  clock. 
At  this  moment  the  subroutine  SAMPLE  is  connected  and  samples  the  analog 
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computer,  storing  the  value  into  ARRAY;  at  the  same  time  it  decreases  a 
counter,  MAC,  which  has  the  effect  of  opening  the  loop  in  SEGSAMP  and 
allows  copying  ARRAY  into  the  proper  location. 

The  sample  is  taken  at  the  circuit  node  connected  to  the  trunk  line 
determined  by  LINE  whose  initial  value  is  one,  but  can  also  be  adjusted 
to  any  legal  trunk  line  designation. 

Figure  2.3  is  a  simplified  flow  chart  that  illustrates  the  operation 
of  the  subroutines. 

D.  SUBROUTINES  DISPLAY  AND  TISPLAY 

These  are  also  the  same  subroutine  with  different  calling  names  to 
allow  their  location  in  two  levels  of  segmentation. 

DISPLAY  is  designed  to  encode  the  functions  as  required  for  display 
on  the  graphics  terminal.  It  has  basically  three  options.  The  first 
option  was  developed  for  the  presentation  of  correlation  and  convolution 
operations.  In  this  case  the  y-axis  is  positioned  at  the  left  of  the 
screen.  The  signals  to  be  correlated  and  the  result  of  the  operation  are 
plotted  in  three  sets  of  coordinate  axes  one  below  the  other.  As  a 
variation  of  this  option  only  the  first  function  and  the  three  coordinate 
axes  or  only  the  other  two  functions  are  displayed  at  one  time;  this 
feature  is  useful  to  the  subroutine  FASTCOR  as  is  shown  later.  The  number 
of  points  used  to  describe  each  function  may  also  be  varied.  The  second 
option  is  applied  to  the  Fourier  transform  and  here  the  same  number  of 
points  are  required  for  all  three  functions.  The  y-axis  is  placed  at  the 
center  of  the  screen.  In  this  case  the  first  one  is  the  time  function 
itself,  the  second  the  absolute  value  of  the  transform,  and  the  third  the 
phase  angle.  The  last  option  is  provided  for  single  function  display. 
The  function  is  presented  at  the  center  of  the  screen  with  the  ordinate 
expanded, 
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FIGURE  2.3 
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E.  SUBROUTINE  FASTCOR 

This  subroutine  constitutes  the  basic  feature  of  the  system  and  its 
purpose  is  to  perform  a  fast  convolution.  Next  will  be  presented  a  review 
of  the  mathematical  concepts  that  were  considered  to  build  this  program. 

The  convolution  of  two  functions  is  defined  as  [Ref.  4  and  8]: 

00  oo 

f(x)  =  r(t)  *  s(t)  =  /  r(x)s(x-x)dx  =  /  r(x-x)s(x)dx 

^oo  —oo 

The  cross-correlation  is: 

00  CO 

R(t)  =  /  rCx)s(x-x)dx  =   /  r(x-x)s(x)dx 

-oo  -co 

This  is  a  similar  integral  to  the  convolution  with  the  difference 

that  now  none  of  the  functions  is  time  reversed. 

If  one  of  the  functions  is  taken,  say  r(t),  and  time  reversed  before 
calculating  its  convolution 

CO 

g(x)  =  rOt)  *  s(t)  =  /  r(x-x)s(x)dx 

-co 

then 

g(x)  =  R(x) 

and  it  is  concluded  that  the  evaluation  of  both  integrals  reduces  to  a 
single  problem,  since  one  of  the  functions  can  always  be  reversed  prior 
to  taking  the  product  and  integrating. 

The  solution  of  these  integrals  is  by  no  means  a  simple  task.  Even 
in  the  case  where  both  functions  are  limited  in  time  so  that  the  integral 
vanishes  outside  a  certain  interval,  that  is,  when  it  is  only  necessary 
to  evaluate 
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x2 

/        r(x)s(x-x)dx 

xi 


the  computation  is  not  straight-forward. 

In  practical   applications  this  integral   is  normally  evaluated  by 
graphical  or  numerical  methods.     In  both  cases  the  approximate  solution 
is  still   time  consuming  in  general. 

Another  method  to  evaluate  the  convolution  integral   can  be  reached 
through  the  theorems  of  Fourier  transformation. 

Fourier  transform  and  its  inverse  are  defined  as: 


F(u>)  =  F[f(t)]  =      /     f(x)e~ja)Xdx 


00 


f(t)  ■  F_1[F(a))]  =^_    /     FMejaJtdo 


-oo 


The  inverse  transform  of  the  product  of  the  transforms  of  each  of 
the  functions  r(t)  ans  s(t)   is  as  follows: 


F~'[RMS(w)]   =     ^_    /     R(aJ)S(a))eJa3Tdco 


00  00 

L    /      [    /     r(x)e"jwxdx]  SCa3)eja3Tda) 


OO  _  OO  •        /  \  °° 

=    /       r(x}dx  [ip    /     SCco)eJajCT'x)da)]  =    /     r(x)s(x-x)dx 

mCO  -CO  -CO 

Since,  as  is  easily  seen,  the  Fourier  transform  of  r(-t)  is  R*(u>)» 

it  immediately  follows  that 

1 
F  [R*(oo)S(w)]  =   /  r(x-T)s(x)dx 
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What  these  expressions  denote  is  that  if  the  Fourier  transform  of 
both  functions  is  known,  the  evaluation  of  the  convolution  or  correlation 
integrals  may  be  obtained  by  multiplying  the  transforms  and  then  performing 
the  inverse  transformation.  This  process  greatly  simplifies  the  task  if 
there  is  a  method  to  easily  compute  Fourier  transforms. 

Such  methods  are  not  known,  in  general.  However  if  the  functions  are 
represented  in  the  discrete  domain,  there  are  fast  solutions  that  can  be 
obtained  with  digital  machines.  The  computation  then  gives  yery   close 
estimates  for  correlation  and  convolution. 

In  the  discrete  case  —  the  functions  are  represented  as  r'(t)  and 
s'(t)  thus  [Ref.  5]: 

00  oo 

r*(t)  =  S  r(nT)   and  s*(t)  =  Z  s(nT) 

n  =  -oj  — oo 

where  T  is  the  spacing  betv/een  samples,  and  if  the  values  of  r(t)  and 

s(t)  are  zero  outside  a  certain  interval,  say  0  to  (N-l)T  as  always 

happens  in  practical  situations, 

N-l  N-l 

r'U)  =  I       r(nT)        s'(t)  =  I       s(nT) 
n=0  n=0 

The  discrete  Fourier  transform  is  defined  as 

N-l 
R(ktt)  =  t       r(nT)e"jfiTnk 
n=0 

with     Q.  =  Tjy    and  an  inverse  transformation 

i     N>1  .nT  . 

r(nT)  =     ^r      z       R(kft)eJ"Tnl< 
N     k=0 
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The  discrete  Fourier  transform  has  several   properties,  many  of  them 
paralleling  those  of  the  continuous  transform.     One  fundamental  distinc- 
tion, however,  results  from  the  definition  of  Q, :     the  transform  or  the 
inverse  transform  has  only  N  distinct  values  for  k(or  n)  between  0  and 
N-l  with  these  values  repeating  themselves  for  k  outside  this  range. 
This  may  be  seen  as  the  transformation  giving  the  result  of  a  periodic 
sequence  with  period  N. 

A  similar  relation  exists  between  the  convolution  and  the  product 
of  the  transform  of  two  functions  as  occurred  in  the  continuous  case: 


N-l 
f(nT)  =  i      Z     R(kfi)S(kn)ejnkTfi 

1N  k=0 

*  1  V   [  V     r(£T)e-J*kTO]     [V     s(mT)e"J'mkTfi]eJ'nkTfi 

IN  k=0      1=0  m=0 

4  V     V     r(£T)s(mT)     [V     eJ(n-nWl)k^ 

"  1=0    m=0  k=0 


It  is  not  difficult  to  show  that  the  summation  inside  the  brackets 
is  zero  unless  m  is  equal  to  (n-£ )  modulo  N  when  the  result  is  zero. 
Hence,  substituting  (n-£)  modulo  N  for  m: 

N-l 
f(nT)  =  I       r(U)s((n-£)T) 
A=0 

where  (n-£)  is  to  be  considered  as  (n-£)  modulo  N.  f(nT)  has  then  an 
expression  similar  to  that  given  for  f(x)  and  it  may  be  said  that  f(nT) 
is  a  convolution.  However  it  has  the  special  characteristic  of  repeating 
itself  with  period  N  and  so  is  named  a  circular  or  periodic  convolution. 
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As  a  conclusion,  an  estimate  of  the  convolution  of  the  functions 

N-l 
r(t)  and  s[t)  is  f'(t)  =  I     f(nT),  where  r  and  s  are  assumed  periodic 

n=0 
with  period  NT. 

To  evaluate  that  expression  it  is  possible  to  take  the  discrete 
Fourier  transform  of  each  function,  multiply  them  together  and  then 
compute  the  inverse  transform. 

Fortunately,  a  s/ery   convenient  and  fast  algorithm  known  as  Fast- 
Fourier-Transform  has  been  developed  and  is  specially  adapted  for  use 
with  digital  machines  [Ref.  3]. 

Thus,  it  is  now  possible  to  evaluate  convolutions  and  correlations 
for  practical,  real-time  applications. 

But  not  all  functions  dealt  with  are  periodic.  On  the  contrary, 
very  often  they  are  limited  to  a  relatively  short  interval  of  time.  The 
method  can  still  be  applied  if  the  functions  are  prepared  before  per- 
forming the  convolution.  Figure  2.4  illustrates  how  this  can  be  done. 
The  convolution  of  two  rectangular  pulses  is  shown  in  Fig.  2.4a  and  in 
2.4b  the  same  thing  is  presented  in  the  case  of  periodic  repetition  of 
the  pulses.  The  periods  were  purposely  chosen  to  space  the  pulse  widely 
so  that,  when  one  series  is  slid  across  the  other,  there  is  no  possibility 
of  any  single  pulse  of  the  first  function  coinciding  with  more  than  one 
pulse  of  the  second  function  at  the  same  time.  One  period  of  the  convo- 
lution has  then  the  same  form  as  the  convolution  of  the  aperiodic  functions 
Therefore,  in  order  to  evaluate  the  aperiodic  convolution  it  is  necessary 
only  to  extend  both  functions  with  zeros,  constructing  periodic  functions 
with  a  period  size  of  M+N,  where  M  and  N  are  the  actual  sample  sizes  of 
the  original  functions,  and  taking  only  one  period  from  the  resulting 
periodic  convolution. 
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Another  application  of  this  method  refers  to  the  case  when  one  of 
the  functions  is  very  short  in  comparison  with  the  other,  as  appears  in 
a  detection  problem  when  a  short  replica  is  correlated  against  a  very 
long  incoming  signal.  This  situation  is  presented  in  Fig.  2.5. 

The  long  sequence  is  partitioned  in  sections  of  size  M,  each  section 
extended  with  N  points  from  the  previous  section,  and  the  samll  one 
extended  with  M-N  zeros.  If  the  correlation  of  the  N  samples  (plus  M 
zeros)  against  successive  N+M  sized  sections  is  performed,  and  the  results 
added,  it  is  evident  that  the  first  N  points  of  each  partial  result  will 
be  originated  in  the  same  way  as  the  last  N  points  of  the  previous  sec- 
tion. This  distortion  is  avoided  if  the  N  first  points  of  e^ery   partial 
result  are  discarded. 

In  the  derivations  mentioned  above,  no  restriction  was  imposed  upon 
the  functions  that  can  in  general  be  complex.  Therefore  in  practice  a 
group  of  four  convolutions  or  correlations  is  performed  at  the  same 
time,  i  .e. , 

r(t)  =  rr(t)  +  jrt(t);  s(t)  =  sr(t)  +  js.(t) 

r(t)*s(t)  =  r  *s  -  r.*s.  +  jr  *s.  +  jr.*s 

In  most  applications,  however,  real  functions  are  dealt  with  and  it 
is  then  possible  to  perform  more  than  one  operation  simultaneously.  For 
instance,  the  partial  result  can  be  evaluated  for  two  sections  in  one 
operation,  or: 

r(t)  =  rr  +  jO    s^t)  =  slr  +  jO    s^t)  =  s2f  +  jO 

ss  ■  slr  +  js2r  r(t)*s$(t)  =  rr*slr  +  jrr*s2r 

After  the  operation  is  concluded  the  real   part  is  taken  as  one  section 
of  the  result  and  the  imaginary  part  as  the  other. 
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The  subroutine  FASTCOR  uses  for  the  evaluation  of  the  Fourier  trans- 
form a  very  fast  program,  F0UR2,  based  on  the  Cooley-Tuckey  method,  which 
is  available  in  the  library  of  the  Computer  Laboratory.  This  program 
requires  that  the  functions  to  be  transformed  have  a  number  of  points 
equal  to  a  power  of  2,  and  is  the  reason  for  this  restriction  on  the 
sample  size. 

FASTCOR  has  three  main  options  either  for  convolution  or  correlation. 
They  are  assigned  by  the  argument  IT,  i.e.,  periodic,  aperiodic,  or 
recursive.  Figure  2.6  shows  a  flowchart  for  FASTCOR. 

Initially  the  size  of  the  samples  is  checked  against  the  option 
selected.  In  any  case,  the  prepared  sample  (a  "period")  for  both  samples 
has  to  have  the  same  number  of  points. 

If  a  periodic  operation  is  chosen  and  the  second  function  is  longer 
than  the  first  one,  with  a  maximum  of  512,  it  is  reduced  to  the  shorter 
length.  If  an  aperiodic  operation  was  selected  and  the  sizes  of  both 
samples  are  equal,  in  which  case  the  maximum  length  is  256  (M) ,  the  two 
functions  are  padded  with  zeros  to  2M.  In  case  the  first  sample  is 
smaller,  it  is  limited  to  one-half  of  its  size  and  the  second  one  has 
the  number  of  points  reduced  to  a  value  equal  to  its  original  value 
minus  one-half  of  the  first  function's  original  size.  Both  functions 
are  then  padded  with  zeros  up  to  the  initial  size  of  the  second  function; 
this  assures  no  overlapping  in  the  result.  Finally,  with  the  recursive 
option  the  first  function  is  always  assumed  to  be  smaller  than  the  second 
which  has  a  maximum  acceptable  length  of  512;  the  first  sample  is  padded 
with  zeros  to  the  length  of  the  second. 

The  next  step  is  to  flip  the  first  function,  send  it  to  be  displayed 
on  the  screen,  since  it  will  soon  be  erased  from  memory,  and  store  it  as 
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complex  with  imaginary  part  zero.  If  a  correlation  is  to  be  performed, 
the  function  is  displayed  without  reversal  but  reversed  and  stored  as 
complex;  for  convolution  the  opposite  occurs. 

If  the  option  is  not  recursive,  the  second  function  is  also  stored 
as  complex  with  zero  imaginary  part,  and  both  are  then  Fourier  trans- 
formed, multiplied  together,  and  the  result  inverse  transformed.  During 
these  operations  some  transfers  of  storage  occur  to  minimize  the  space 
used.  Finally,  the  second  function  and  the  result  are  shown  on  the 
screen  unless  no  display  is  required. 

When  the  option  is  recursive,  an  auto-correlation  is  initially  per- 
formed on  the  first  function  to  establish  the  normalization  coefficient 
and  the  actual  threshold  value.  This  normalization  factor  is  computed, 
since  the  maximum  amplitude  a  correlation  can  attain  occurs  when  the 
two  functions  are  identical  and  therefore  corresponds  to  the  maximum 
value  of  the  auto-correlation.  But  since  the  amplitudes  of  the  two 
functions  have  been  previously  and  independently  normalized,  the  new 
normalization  factor  is  corrected  by  the  ratio  Q  mentioned  earlier.  It 
was  chosen  to  normalize  the  threshold  value  to  the  maximum  amplitude  of 
the  auto-correlation,  which  means  that  a  relative  threshold  of  one  would 
correspond  to  a  perfect  matching  of  the  functions.  This  actual  value  of 
threshold  is  then  further  normalized,  in  the  same  way  as  the  result,  for 
display  purposes. 

Having  completed  this  task,  the  second  function  is  then  stored  as 
complex.  In  the  strictly  recursive  case  a  double-sized  (2M)  sample  is 
stored  as  real,  and  the  next  section,  with  overlap,  as  imaginary.  For 
the  determination  of  probability  of  detection,  however,  a  sample  section 
padded  with  zeros  is  stored  as  real ,  and  the  next  section  as  imaginary 
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so  that  effectively  a  double  aperiodic  correlation  is  taken  at  one 
time.  After  correlation,  the  two  sections  are  moved  in  core,  compared 
against  the  threshold,  and  displayed  if  so  desired. 

If  a  correlation  is  being  evaluated,  the  threshold  level  is  also 
displayed,  and  an  indication  appears  as  to  whether  this  amplitude  is 
reached  or  not.   Information  also  comes  to  the  screen  about  the  probabil- 
ity of  detection  that  has  meaning  only  if  this  option  is  in  execution. 
The  last  measured  input  signal-to-noise  ratio  is  also  shown. 

Every  time  the  threshold  is  met  in  a  distance  measuring  operation, 
the  actual  distance,  represented  by  the  sample  number,  is  displayed. 
The  signal,  in  this  case,  is  added  to  noise  only  at  one  predetermined 
distance,  i.e.,  only  in  one  group  of  samples. 

After  this  phase  is  completed  the  first  2M  samples  of  the  second 
function  (signal)  are  discarded  and  the  last  M  are  translated  to  consti- 
tute the  new  first  set.  SEGSAMP  is  then  called  with  INDEX  equal  to  3 
to  fill  out  the  function  with  2M  new  samples,  and  the  computation  starts 
again. 

Exit  from  the  recursion  loop  results  when  either  the  loop  has  been 
executed  a  preset  number  of  times  or  when  the  operator  presses  a  switch 
at  the  graphics  console. 

In  all  cases  the  control  returns  automatically  to  the  sequence 
indicated  by  LABEL  except  when  the  probability  of  detection  (or  the  false 
alarm  rate)  was  requested  in  the  ANALYSIS  mode.  In  that  case  the  proba- 
bility of  detection  (or  the  false  alarm  rate),  threshold  setting,  number 
of  sections  processed  and  input  signal-to-noise  ratio  are  all  printed 
on  the  line  printer,  and  the  operator  must  command  continuation  by 
pressing  the  RETURN  key  of  the  terminal  keyboard. 
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F.  SUBROUTINES  ANALYSIS  AND  SEGFOUR 

These  subroutines  belong  to  the  ANALYSIS  mode. 

ANALYSIS  is  intended  to  clarify  further  the  request  of  the  operator 
and  route  the  computer  to  the  appropriate  sequence.  The  argument  IENT 
is  a  pointer  to  indicate  whether  the  subroutine  is  entered  to  initiate 
a  sequence  or  to  end  it.  This  is  a  requirement  of  the  overlay  structure 
since  the  space  it  occupies  in  core  has  to  be  released  whenever  FASTCOR 
is  called. 

The  initial  selection  is  indicated  by  IJ,  kept  in  common  storage. 

IJ  equal  to  one  indicates  that  sampling  is  desired.  Before  calling 
SEGSAMP  a  reminder  is  displayed  asking  the  operator  if  the  constants 
are  adjusted  to  sample  the  correct  function,  for  the  correct  processor 
and  from  the  correct  point  in  the  analog  (INDEX,  ISAMP,  LINE).  If  the 
answer  is  negative  the  constants  are  displayed  to  allow  corrections 
after  which  sample  has  to  be  requested  again.  An  affirmative  answer 
initiates  the  sampling  operation  and  then  the  new  sample  is  saved  in 
background  storage. 

IJ  equal  to  two  indicates  Fast-Fourier-Transform  and  the  operator 
is  requested  first  to  indicate  which  function:  replica  or  M  sized 
function,  signal  or  2M  sized  function,  correlation  or  the  result  of  the 
last  operation  performed  by  FASTCOR.  Then  SEGFOUR  is  called  to  prepare 
the  function  selected,  perform  the  transformation  through  F0UR2,  calcu- 
late the  absolute  value  and  the  phase  angle  of  the  transform,  send  them 
and  the  original  function  to  the  display  and  read  back  from  secondary 
storage  the  functions  destroyed  during  the  operation. 

When  correlation  is  selected,  IJ=3,  it  is  only  necessary  to  indicate 
if  the  mode  desired  is  periodic,  aperiodic  or  recursive.  FASTCOR  is  now 
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called  through  FASTBRG,  a  bridge  subroutine  required  to  save  addresses 
and  parameters  while  the  segment  swapping  takes  place.  Back  from  FASTCOR, 
FASTBRG  re-enters  ANALYSIS  to  restore  the  functions  destroyed  and  put  the 
result  of  the  correlation  in  backstorage. 

For  IJ  equal  to  four  (convolution),  five  (power  spectrum),  or  seven 
(auto-correlation),  the  operator  must  select  the  function  and  mode.  With 
convolution,  the  "other  function"  is  always  assumed  to  be  the  replica; 
with  the  other  two  operations,  recursive  is  an  illegal  choice  and  if 
made  results  in  an  error  message.  The  functions  are  then  moved  in 
storage  so  as  to  send  to  FASTCOR  the  appropriate  parameters.  Returning 
to  ANALYSIS,  action  similar  to  that  for  correlation  is  taken,  except  when 
power  spectrum  is  the  operation. 

The  evaluation  of  power  spectrum  is  done  through  the  application  of 
the  Wiener-Kinchin  relations  [Ref.  8]. 

It  was  earlier  recalled  that  the  convolution  is  the  inverse  Fourier 
transform  of  the  product  of  the  transform  of  the  functions  to  be  convolved 
and  that  the  correlation  can  be  treated  as  the  convolution  taken  with 
one  of  the  functions  time-reversed  prior  to  the  operation.  Also  the 
transform  of  r(-t)  is  equal  to  the  complex  conjugate  of  R(oj).  Now, 
the  auto-correlation  of  a  signal  for  which  the  energy  content  is  not 
finite  is  defined  as 

T 

pM  =  fj£  \    f  2  r(t)r(t-x)dt   . 

T 


Hence, 


PCt)  =  f&  |[r(t)*r(-t)] 
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Jlim 


PCt)     =     T^    WT       f     RCo3)R*(a))eJWTda) 


f£      /       |RC2Trf)|2ej27TfTdf 


—  CO 


gg 

=     /     G(f)  ej27TfTdf 


where  G(f)  is  the  power  spectral  density.  If  the  energy  content  is 
finite,  the  energy  spectrum  is  considered  instead  of  the  power  spectral 
density.  Therefore,  since  the  power  (or  energy)  spectrum  and  the  auto- 
correlation form  a  transform  pair,  the  spectrum  may  be  evaluated  by 
computing  the  auto-correlation  and  taking  its  Fourier  transform. 

ANALYSIS  then  determines  the  auto-correlation  to  be  performed  and 
when  re-entered,  calls  SEGFOUR  to  transform  the  result.  (None  of  these 
results  are  saved. ) 

The  last  control  of  the  subroutine  corresponds  to  IJ  equal  to  six 
or  the  display  selection.  After  the  answer  for  the  function  is  selected 
and  decoded  the  subroutine  DISPLAY  is  called  to  show  that  function. 

ANALYSIS  has  to  keep  track  of  the  number  of  samples  related  to  each 
selection  and  controls  the  allowed  limit  for  all  situations  reducing  the 
size  of  the  sample  if  necessary.  The  common  location  II  is  used  for 
this  purpose  during  the  overlays  substitution.  The  size  of  the  last 
result  (CORR)  stored  is  saved  in  I MEM.   Its  value  is  set  when  the  sub- 
routine is  re-entered  and  corresponds  to  the  number  of  points  displayed 
as  indicated  by  LBG(3). 

The  flowchart  for  this  subroutine  is  shown  in  Fig.  2.7. 
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G.  SUBROUTINES  PROCESS  AND  PISPLAY 

PROCESS  is  constructed  to  chain  various  sequences  in  order  to  execute 
the  operations  necessary  to  compute  input  signal-to-noise  ratio,  evaluate 
the  processor  input-output  performance,  and  determine  the  ROC  curves.   It 
also  permits  simulating  a  situation  of  target  detection  and  distance 
measurement.  This  automatic  procedure  constitutes  a  part  of  the  SIMULATION 
mode  except  for  the  signal-to-noise  calculation  that  may  be  called 
separately  in  the  ANALYSIS  mode. 

Except  for  distance  measurement,  the  program  establishes  loops  where 
some  parameters  are  varied  for  ey^ry   recursion  until  a  limit  is  reached. 
Programwise,  however,  these  loops  are  interrupted  several  times  by  trans- 
fer to  another  overlay.  More  common  storage  is  then  required  to  save  the 
status  during  the  segment  swapping,  and  to  allow  resumption  of  the  loop 
in  the  proper  position. 

A  simplified  flowchart  for  these  sequences  is  presented  in  Fig.  2.8. 

When  signal-to-noise  ratio  is  called  in  the  ANALYSIS  mode,  it  is 
initially  necessary  to  scan  the  signal  and  noise  injection  potentiometers, 
save  their  settings,  and  restore  them  after  the  computation  is  completed. 
Entered  the  first  time  for  performance  or  ROC  determination,  it  determines 
the  number  of  steps  that  will  be  used  to  calculate  the  SNR.  The  operation 
really  starts  by  saving  the  processor  status  and  establishing  a  linear 
processor,  since  the  input  signal-to-noise  ratio  (after  the  transducer) 
is  desired;  the  processor  is  reset  after  the  computation.  Two  loop  paths 
are  then  executed,  one  to  determine  the  signal  power  at  that  level  of 
injection  (no  noise  present),  and  the  other  for  noise  power  (no  signal 
present).  The  ratio,  in  decibels,  is  taken  and  stored  as  SNR,  and 
displayed  if  the  operation  CONST  is  selected  next. 
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The  Wiener-Kinchin  relation  again 


P(t)     =     /     G(f)  ej27TfTdf 


When  x  is  zero 

00 

P(0)     =     /     G(f)  df 

— oo 

meaning  that  the  value  of  the  autocorrelation  when  x  =  0  is  the  total 
average  power. 

The  absolute  value  of  the  autocorrelation  is 

|P(t)|     ■       |/     G(f)ej2TrfTdf|     <     /       G(f)eje7TfTdf|     =     /     G(f)df     =  P(0). 

—00  .00  —00 

So,  the  power  content  can  be  found  by  taking  the  maximum  value  of  the 
autocorrelation  since  P(0)  is  an  absolute  maximum.  This  is  exactly  what 
the  program  does . 

The  input-output  performance  or  processor  gain  can  be  defined  by 
the  equation 

SNR  =  f(p,W)SNR. 

where  SNR  is  peak  output  signal  power  to  output  noise  power  ratio 
SNR.  is  the  input  signal-to-noise  power  ratio,  and  f(p,W)  is  a  function 
of  the  processor  and  its  associated  information  bandwidth. 

The  performance  operation  begins  sampling  and  storing  the  replica 
in  a  position  never  to  be  destroyed  during  the  process  in  this  particular 
case.  Then  it  uses  the  previously  established  values  for  signal  and 
noise  injection  that  correspond  to  the  highest  input  signal-to-noise 
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ratio  to  be  considered,  and  sequentially  computes  input  signal-to-noise 
ratio,  samples  noise  alone,  cross-correlates   replica  and  noise,  auto- 
correlates  this  output  and  takes  the  maximum  (output  noise  power).     It 
then  samples  signal   plus  noise,  cross-correlates  replica  and  signal   plus 
noise,  finds   the  maximum  (signal   peak  output).     Next  it  reduces  the 
signal   injection  and  correspondingly  increases  the  noise  by  a  pre- 
determined step  and  recycles  as  long  as  the  signal-to-noise  ratio  is 
greater  than  zero.     During  the  operation  the  individual   values  of  signal- 
to-noise  ratio,  output  noise  power,  and  signal   peak-output  are  placed  in 
common  storage.     When  the  loop  ends,  the  computed  values  are  rearranged 
in  core  and  SNR     is  computed.     Finally,  the  subroutine  PISPLAY  is  called 
to  display  the  resulting  curve  with  the  abscissa  as  SNR. (db)  and  the 
ordinate  as  SNR  (db).     (The  input  signal-to-noise  ratio  is  measured  after 
the  transducer  filter.) 

One  detail  must  be  examined  in  the  output  noise  power  computation. 
Since  the  cross-correlation  of  replica  and  noise  has  2M  points,  it  may 
exceed  the  maximum  limit  allowable  for  autocorrelation  if  M  is  256.     The 
solution  adopted  was  to  calculate  the  autocorrelation  in  two  steps, 
one-half  at  a  time  and  then  add  the  two  maxima;  this  is  legitimate  since 
the  correlation  is  a  linear  operation. 

For  ROC  determination,  the  procedure,  up  to  the  input  signal-to-noise 
ratio  evaluation  and  the  recursion  for  the  several   levels  of  signal,   is 
identical   to  the  sequence  just  described,  with  the  exception  that  the 
replica  is  placed  in  secondary  storage.     Inside  the  loop  however  two  other 
loops  exist:     the  outer  one  for  the  threshold  value  that  steps  down  in 
0,1   steps  from  the  initial   setting  until   zero  is  reached;  the  innermost 
loop  has  two  paths,  one  for  the  calculation  of  the  probability  of 
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detection  (signal  plus  noise)  and  the  other  for  false  alarm  rate  (noise 
only).  PISPLAY  is  also  used  to  present  the  results  on  the  screen.  The 
abscissa  is  presented  as  false  alarm  probability  and  the  ordinate  is 
false  alarm  rate  (per  sample).  A  listing  of  the  signal-to-noise  ratio 
for  the  curves  presented  and  normalized  threshold  are  also  displayed. 

According  to  this  description,  instead  of  the  receiver  operating 
characteristic  the  result  is  a  modified  ROC  where  false  alarm  rate  was 
substituted  for  false  alarm  probability.  This  modification  has  more 
practical  significance  and  is  also  easy  to  evaluate  in  a  short  period 
of  time  while  the  probability  of  false  alarm  is  not.  These  two  param- 
eters are  related  by  the  noise  output  bandwidth. 

At  this  point  a  warning  should  be  made.  The  above  estimations  serve 
as  an  indication  of  the  probable  behavior  of  a  processor  and  that  is  the 
purpose  of  the  system.  More  precise  ensemble  average  estimations  would 
require  lengthy  repetitions,  and  it  was  not  considered  worthwhile  to 
incorporate  them.  However,  as  will  be  pointed  out  later,  the  basic  idea 
can  be  used  for  a  much  less  flexible  fast  analytical  system  or  for 
special  purpose  designs. 

Detection  and  distance  measurement  form  the  last  option  of  PROCESS 
which  sets  parameters  and  lines  to  initiate  the  procedure,  and  samples 
the  replica  and  the  initial  sections  of  signal.  A  line  is  set  to  dis- 
connect the  pattern  generator.  I COUNT  is  set  to  allow  a  maximum  distance 
scan  corresponding  to  20M  samples  and  ISTF  to  connect  the  signal  at  a 
distance  of  12M  sampling  intervals. 

H.  SUBROUTINE  CIRCPAT 

The  SIMULATION  mode  is  completed  with  this  subroutine  which  is  able 
to  respond  to  the  typed  commands  and  set  circuits  and  parameters  to 
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simulate  different  transmission  signals,  such  as  CW  burst,  FM  slide, 
FSK  and  propeller  signature. 

These  simulations  are  better  discussed  in  the  section  dealing  with 
applications. 

I.  AUXILIARY  SUBROUTINES 

A  number  of  auxiliary  subroutines  were  developed  to  enhance  the 
utility  and  convenience  of  this  system  for  the  investigation. 

GRAPHAL  has  already  been  referred  to  and  permits  the  control  of  the 
analog  computer  from  the  graphics  terminal.  It  is  entered  whenever  the 
command  for  one  of  its  functions  is  typed;  if  this  function  is  to  scan 
the  analog,  ANOSCAN  is  called.  These  subroutines  may  be  modified  for 
general  use. 

PPLOT  copies  the  functions  displayed  on  the  screen  and  has  three 
options  to  output  them.  The  first  is  through  the  library  routine  VPLOT, 
with  identifier  changed  to  WPLOT  since  it  was  necessary  to  have  it  loaded 
only  in  the  overlay  in  order  to  reduce  core  storage  requirements.  It 
plots  the  curves  on  the  line  printer  in  a  one-page  sized  plot  where  the 
x-axis  is  transverse  to  the  paper.  NONGPLOT,  the  new  identifier  given 
to  the  library  subroutine  LONGPLOT  for  the  same  reason  as  above,  also 
prints  on  the  line  printer  but  with  the  x-axis  longitudinal  to  the  paper 
and  a  plot  occupies  a  variable  number  of  pages  depending  on  the  number 
of  points  plotted  and  on  a  parameter  set  by  the  operator.  Through  another 
sequence,  PPLOT  uses  the  eight  channel  recorder  as  the  output  medium. 
When  selected,  it  first  asks  that  the  recorder  be  made  ready,  and  when  the 
operator  presses  the  RETURN  key  it  waits  20  seconds  before  starting  the 
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plot.  The  information  about  the  potentiometer  P0I1  displayed  there 
indicates  the  normalized  threshold  level  that,  if  desired,  may  be  intro- 
duced in  that  potentiometer.  For  this  case,  an  analog  probability  counter 
will  deliver  pulses  to  the  recorder  whenever  the  correlation  curve  crosses 
the  threshold  value.  The  value  for  adjustment  of  P 0 11 ,  however,  is 
meaningful  only  if  the  curve  displayed  corresponds  to  a  recursive  type 
of  correlation. 

Subroutine  DEMANAL  routes  the  command  to  other  subroutines  which 
perform  correlations  through  a  direct  approach,  that  is,  by  simple  inter- 
pretation of  the  integral  definition.  Those  subroutines  are:  ANACOR 
which  performs  a  periodic  autocorrelation  using  an  analog  multiplier  and 
integrator;  DIGCOR  which  evaluates  auto  or  cross-correlation  entirely 
digitally;  and  RECOR  which  computes  recursive  cross-correlation. 

Subroutine  NORMALIZE  is  intended  to  normalize  a  function  to  some 
preset  value.  The  storage  address,  the  maximum  value  desired  and  the 
number  of  points  are  the  arguments  passed  to  it. 

When  the  program  is  first  loaded  or  after  GRAPHAL  is  called,  INITIAL 
initializes  the  text  and  graph  directory  and  the  graphics  terminal.  It 
is  entered  to  erase  the  text  blocks  after  an  operation  is  concluded  or 
to  erase  the  graph  blocks  after  a  new  selection  is  made  (except  if  the 
selection  is  CONSTANTS  or  any  of  the  plot  options). 

J,   SUBROUTINES  DECO  AND  CONST 

DECO  works  as  the  control  section  of  the  system,  allowing  the  user 
to  type  in  his  selection  which  is  then  decoded  and  interpreted  to  channel 
the  actions  in  the  desired  direction.  It  forms  a  table  of  control  words, 
composed  of  the  first  four  letters  of  the  selections  available,  and  when 
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a  command  is  decoded  checks  it  against  that  table  and  sets  two  pointers, 
II  and  IJ,  and  returns  to  MAIN.  II  is  then  used  to  select  the  segment 
while  IJ  indicates  the  correct  choice.  If  no  entry  exists  in  the  table 
matching  the  selection,  an  error  message  is  displayed  and  a  new  order 
may  be  entered.  DECO  is  called  automatically  after  the  completion  of 
an  operation  or  when  an  illegal  procedure  is  requested  in  another  sub- 
routine. 

The  control  function  is  completed  with  CONST  which  displays  various 
parameters  and  pointers  and  permits  the  operator  to  alter  them  inside 
allowable  limits.  The  names  under  which  they  appear  are  not  necessarily 
those  used  in  the  program,  but  others  more  easily  remembered.  A  full 
discussion  of  this  will  be  given  with  the  graphics  terminal  usage 
explanation. 
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III.  ANALOG  CIRCUITS 

The  analog  computer  was  selected  to  simulate  the  signals  because 
they  are  generated  in  the  same  form  as  they  occur  in  reality.  Thus  the 
analysis  of  existing  signals  can  be  made  by  the  simple  substitution  of 
the  actual  signals.  In  this  case  the  analog  computer  itself  acts  as 
the  interface  to  the  system. 

Another  important  feature  considered  was  the  relatively  large 
parallel  logic  capability  of  the  COMCOR  Ci-5000  analog  computer,  which 
allowed  construction  of  control  circuits  to  enlarge  the  range  of  utili- 
zation and  analog  processing. 

The  analog  computer  circuits  included:  one  sine  generator  used  as 
PATTERN  GENERATOR,  one  HARMONIC  GENERATOR  fed  by  three  independent  sine 
generators,  one  WHITE-NOISE  generator,  a  SIGNAL  COMBINER,  a  BANDPASS 
FILTER,  an  INTERRUPT  GENERATOR  and  some  auxiliary  circuits. 

As  a  general  rule  the  circuits  and  their  parameters  are  remotely 
controlled  from  the  graphics  terminal  with  exception  of  those  few  which 
need  patch  relocation.  All  lines  that  can  be  set  are  tied  to  another 
test  line  in  order  to  permit  status  check  from  the  display. 

A.  PATTERN  GENERATOR 

1 .  Theory 

..   .2     2 
The  differential  equation  y  =  — %■  =  -00  y  represents  a  sine  wave 

dr 

and  its  solution  can  easily  be  implemented  with  integrators  and  adders. 
The  mechanization  of  this  equation  is  presented  in  the  following  expres- 
sions where  the  bar  means  that  the  variable  is  a  scaled  (computer) 
variable,  them's  are  amplitude  scaling  coefficients,  B  is  a  time  scaling 
coefficient,  and  t  is  comoputer  scaled  time  interval: 
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2.  Implementation 

The  basic  implementation  of  the  equations  above  is  shown  in 
Fig.  3.1. 

Since  it  was  desired  to  make  this  circuit  simulate  a  CW  burst, 
as  well  as  a  linear  FM  slide  and  an  FSK  codes,  a  multiplier  was  utilized 
for  p-,  and  p2  because  it  allows  changing  frequency  during  operation 
without  resetting  the  amplifiers.  Additional  circuitry  was  required  to 
generate  the  multiplier  reference  voltage,  to  select  the  type  of  opera- 
tion, and  to  switch  the  generator  on  or  off.  The  final  diagram  is 
presented  in  Fig.  3.2.  A017  is  the  integrator  for  y  and  A021  for  -y. 
A203  and  A204  provide  polarity  reversal  as  required.  The  multiplier 
M401  implements  p,  and  M402  introduces  p„,  Relay  R26  connects  the 
reference  input  of  the  multiplier  either  to  the  ramp  generator  A025  or 
to  the  FSK  relay  RIO.  P030  and  P031  adjust  the  lower  and  upper  shift 
reference  using  A205  and  A206  as  loads.  R006  is  the  relay  that  switches 
the  generator  output  on  or  off. 

The  ramp  generator  computes  the  value  of  R  =  I  +  St  where  R  is 
the  reference  voltage,'  I  the  initial  setting  and  S  the  slope.  P026 
establishes  I,  P025  fixes  S  and  P027  fixes  the  threshold  value  for  C07 
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which  controls  the  reset  time  for  the  ramp.  P032  and  C04  are  responsible 
for  the  shift  instant.  P021  prevents  the  necessity  to  set  P025  to  very 
small  settings  and  A027  and  A030  serve  as  load  to  A025  when  R26  is  open. 

To  make  it  possible  to  run  this  generator  at  very   high  actual 
frequencies  [150  Hz),  positive  feedback  is  introduced  by  P017.  The 
setting  of  this  potentiometer  is  not  critical  since  the  limiter  LOO 
prevents  the  amplitude  from  growing. 

3.  Logic 

The  logic  circuitry  required  to  control  the  pattern  generator, 
presented  in  Fig.  3.3,  is  the  result  of  the  logic  equations  developed 
in  the  Appendix. 

Lines  T007  and  T010  are  set  to  control  the  integrator  mode 
through  several  drivers. 

The  NAND  gates,  comparator  C07,  delay-flop  DF02  and  line  T011 
direct  the  ramp  generator  operation.   It  then  follows  the  general  mode 
except  if  CW  is  required  when  it  remains  always  in  reset.  During  compu- 
tation the  ramp  generator  is  made  to  reset  every  time  a  sample  period 
ends. 

C04  actuates  the  FSK  switch  RIO.  Lines  T016  and  T017  bring  the 
voltage  to  R26  and  R06,  respectively. 

4.  Scaling 

The  presentation  following  has  the  purpose  of  establishing  simple 
reference  tables  for  quick  problem  scaling.  The  final  objective  is  to 
find  the  value  of  R,  the  reference  voltage,  and  the  amplifier  gain  suit- 
able for  each  case. 

The  amplifier  gain  (G)  is  a  function  of  the  input  conductance 
(I-mmhos)  and  the  feedback  capacitor  (C-mfd),  coefficients  p,  and  p^ 
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mentioned  earlier  are  the  product  R6  divided  by  the  constant  100, 
characteristic  of  the  multiplier.  For  convenience,  R  was  multiplied 
by  100  in  the  final  result  to  become  the  integer  value  that  is  sent  to 
set  the  potentiometers. 

n  -  -BJL  -  2L  .   r  -  *c  -  f  -  u 

Hi   innr   ™a  »   p 


1   100C   aS  '      t    f    oo 

c   c 

where  a  and  a  are  amplitude  scaling  factors  for  y  and  y,  B  the  time 

scale  factor,  t  the  computer  scaled  time  interval  and  f  and  u  the 
c       r  c     c 

scaled  frequency  and  angular  frequency  respectively. 

Since  y  =  A  sin  cot,  y  =  -Aco  cos  tot  assuming  A  equal  100  for  a 
maximum  voltage  scaling  a  =  1  and  a  =  to.  Thus 

p  =  100  x  100  C  oo   20000TTf  C 
K       61  31 

R  has  to  be  kept  below  10000  (lOOv)  and  it  was  chosen  to  keep 
C  as  .1  for  real  time  and  .01  for  time  scale  for  an  actual  computer 
frequency,  f  ,  around  15  Hz  and  150  Hz  respectively.  The  other  values 
follow  as  a  consequence  and  are  tabulated  in  Table  I  in  the  Appendix. 
5.  Operation 

For  CW  burst  code,  A027  is  kept  in  RESET  through  T011,  R26  is 
held  in  the  normally  closed  position  (NC) ;  P025  is  set  to  zero  and  P026 
to  the  reference  voltage  R. 

In  FM  slide  A027  is  made  to  compute  once  every  sample  and  then 
reset.  The  relay  R26  is  maintained  in  the  normally  closed  position(NC) 
The  potentiometers  are  set  in  the  following  way: 

P021  to  the  inverse  'of  the  gain  of  A025; 

P025  to  the  ratio  of  the  upper  reference  voltage(Ry)  minus  the  lower 
reference (R[_)  to  the  actual  sampling  time,  i.e., 
(R[j-R|_)  x  sampling  period  x  number  of  points  in  the  sample. 
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P026  to  the  lower  reference  voltage;  and 
P027  to  the  upper  reference  voltage. 

For  FSK  the  same  adjustments  as  for  FM  slide  are  necessary  with 

the  following  alterations: 

R26  -  placed  in  normally  open  (NO)  through  the  action  of  T016; 

P030  -  to  the  lower  reference  voltage; 

P031  -  to  the  upper  reference  voltage; 

P032  -  to  the  central  value  of  the  reference  voltage  R. 

In  all  cases  above,  if  distance  measurement  is  required  the 

relay  R06  is  maintained  in  normally  open  (NO)  all  the  time  except  at  the 

moment  when  the  burst  is  to  occur  when  it  closes  under  control  of  T017. 

B.  HARMONIC  GENERATOR 

1 .  Sine  generators 

The  three  sinusoidal  oscillators  included  in  the  harmonic  generator 
are  of  the  basic  type  presented  in  Fig.  3.1,  with  p,  and  p?  adjusted  by 
potentiometers.  The  scaling  formulas  for  the  multiplier  case  also  apply 
here. 

The  harmonic  generator  was  built  to  simulate  a  relatively  low- 
pass  line  spectra  similar  to  the  signal  generated  by  a  propeller.  One  of 
the  sine  generators  produces  the  signature  while  the  others  produce  a 
"non-white"  noise,  as  if  they  were  other  identical  sources  of  slightly 
different  frequency. 

2.  Multipliers 

In  the  diagram  of  Fig.  3.4  it  is  seen  that  the  three  generators 
were  combined  and  fed  to  four  multipliers  arranged  to  give  as  output, 
after  combination  in  A052: 


E  =  ^e  +  K2e2  +  K3e3  +  K4e4  +  K5e5  +  K6e6 
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The  parameters  K  are  adjusted  in  order  by  P051  to  P056  with  the  first 
three  feeding  the  amplifier  at  I  =  10  and  the  last  at  I  =  1 .  This 
requires  that  P051  through  P053  have  a  maximum  of  10  while  the  sum  of 
all  three,  times  10,  plus  the  remaining  three  be  10000  or  less.  Any 
one  of  the  parameters  may  have  its  polarity  reversed  by  the  switches 
R00,  R02,  R04,  R20,  R22 ,  R24  activated  by  the  lines  T001  through  T006. 

Table  II  in  the  Appendix  gives  some  of  the  combinations  of  coef- 
ficients and  frequencies  up  to  sixth  harmonic  when  one  or  two  frequencies 
constitute  the  input  signal. 

C,  WHITE  NOISE  GENERATOR 
1 .  Theory 

To  build  a  source  of  random  noise  a  pseudo-random  sequence 
generator  was  used. 

As  is  well-known,  a  true  random  binary  sequence  may  be  classified 
by  three  properties  [Ref.2j: 

a)  The  number  of  ones  approximately  equal  to  number  of  zeros; 

b)  About  half  of  the  no-change  states  have  the  length  of  one, 
about  one-fourth  have  length  two,  one-eighth  have  length  three, 
and  in  general  l/2n  have  the  length  n; 

c)  The  expected  value  of  the  autocorrelation  has  a  peak  at  the  origin 
and  decreases  rapidly  away  from  there. 

A  linear  maximal  sequence,  like  those  generated  by  specially 

designed  feedback  shift  registers,  has  properties  yery   close  to  these: 

a)  The  number  of  ones  is  exactly  one  more  than  the  number  of  zeros; 

b)  The  number  of  the  no-change  states  of  length  n  is  l/2n; 

c)  The  finite  autocorrelation,  the  sequence  considered  periodic, 
has  the  form  of 

L 

P(k)  =  Z  a  a  .,  ,  where  L  is  the  sequence  length  and  an  the 
v  '     ,  n  n+k  ^        ^       n 

th 
value  of  the  n   point  of  the  sequence.  This  expression  has  only 

two  values,  say  L  when  k  =  0,  ±L,  ±  2L,  . . .  and  -1  otherwise. 
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Thus,  within  appropriate  limits,  a  pseudo-random  sequence  may 
be  used  to  simulate  a  true  random  one. 

To  generate  binary  analog  pseudo-random  noise  from  a  random 
sequence,  two  techniques  can  be  used.  In  one  of  them  all  bits  are  added 
modulo  2  for  each  shift.  For  the  other  only  one  bit  is  taken  in  each 
shift.  This  logic  signal  is  then  amplified  and  clipped  to  originate  a 
pseudo-random  square  wave.  The  second  technique  was  adopted. 

The  analog  noise  generated  has  first  and  second  order  statistics 
essentially  identical  to  those  of  true  random  noise  if  the  time  of 
observation  is  short  compared  with  the  sequence  period.  Also,  it  is 
found  in  the  literature  that,  after  lowpass  filtering,  high  confidence 
level  tests  indicate  essentially  a  first-order  Gaussian  behavior  [Ref.  7]. 

Again,  for  a  long  sequence  and  a  relatively  short  observation 
interval,  it  is  not  possible  to  resolve  the  spectral  lines  that  neces- 
sarily exist  in  the  power  spectrum  of  a  periodic  waveform. 
2.  Implementation 

As  shown  in  Fig.  3.5,  a  25-bit  feedback  shift  register  (FSR)  was 
constructed  with  bits  3  and  25  being  added  modulo  2  by  the  quarter-adder 

formed  by  A400,  A200,  A100,  and  then  fed  back  to  the  input.  This  arrange- 

25 
ment  generates  a  maximal  sequence  with  length  L  =  2   -  1  =  33,  554,  431. 

The  output  of  the  FSR  and  its  inverse  were  used  to  actuate 
digital -analog  switches,  DA01  and  DA02,  which  apply  -lOOv  or  +100v, 
respectively,  to  the  amplifiers  A010  and  A016,  whose  outputs  are  further 
combined  in  A014. 

In  order  to  avoid  the  all -zero  state,  bit  one  is  set  at  reset 
time  and  all  the  remaining  reset. 

One  requirement  of  a  system  built  for  signal  analysis  is  the 
possibility  to  return  to  any  observed  situation  under  operator  control. 
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For  this  reason,  the  shift  pulse  and  the  reset  lines  were  made  logically 
controlled  by  a  special  circuit,  to  be  described  later,  so  that  the  FSR 
operation  is  synchronized  to  the  generator  modes. 

D.  SIGNAL  COMBINER 

All  three  types  of  signals  just  described  were  brought  to  an 
amplifier,  A024,  through  the  potentiometers,  P022,  P023,  P024,  as  shown 
in  Fig.  3.6.  The  amount  of  injection  of  each  one  can  be  controlled  and 
any  combination  achieved. 

Potentiometer  P015  was  included  to  feed  an  initial  condition  to  A015, 
if  desired.  With  this  feature  and  no  inputs  to  A015  and  A024,  it  is 
possible  to  simulate  an  impulse  to  the  filters  that  follow  and  then 
analyze  their  characteristics. 

Line  T500  was  connected  at  the  output  of  the  combiner  to  make  it 
possible  to  sample  before  the  input  filter  if  desired. 

E.  BANDPASS  FILTER 

Actually,  two  filters  were  included  after  the  combiner  to  simulate 
the  transducer  and  the  input  circuitry.     The  original   intention  was  to 
have  a  narrow-bandpass  filter  for  the  active  detection  simulation  and 
the  lowpass  filter  for  the  propeller  signature  problem. 

1 .     Narrow-band  Filter 

A  second-order  filter,  shown  in  Fig.   3.7  was   installed  for  this 
purpose.     Its  frequency  response  is  given  by 

JojA 


2     2 

(to     -co   )   +  jooB 


to 
where  A  is  the  gain  and  B  the  bandwidth  related  to  the  Q  by  B  s  —  . 
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The  differential  equation  describing  this  circuit  is: 

y  +  ABy  +  Ao>     y  =  -Ax 
which,  when  integrated  and  put  in  terms  of  computer  variables ,  becomes 

y=-/|ydtc-/^zdtc-/ixdtc 


z  =  -ff-  C-y)  dtc 


P036  then  serves  to  adjust  the  gain,  P037  the  bandwidth,  and 
P040  and  P041  the  central  frequency.  The  values  of  Table  I  may  also 
be  used  for  these  potentiometers  if  the  frequency  is  substituted  by 
3A/2tt,  Af  and  fQ  respectively.  To  avoid  small  settings  for  P036  this 
potentiometer  was  normally  connected  to  an  amplifier  having  a  gain  ten 
times  smaller  than  that  used  for  the  f«  input. 

The  filter  characteristic  may  be  estimated  by  application  of  an 
impulse,  as  just  mentioned.  An  example  of  this  was  taken  and  is  shown 
in  Fig.  1.9. 

Line  T501 ,  connected  at  the  filter  output,  makes  it  possible 
to  take  samples  at  this  point. 
2.  Lowpass  Filter 

A  simple  RC  filter  was  used  for  lowpass. 

To  assure  the  first-order  Gaussian  characteristic,  the  time 
constant  RC  has  to  be  much  greater  than  one  shift  period,  At,  but  still 
less  than  25At,  the  sequence  period.  Of  course,  the  inverse  of  RC  is 
the  cutoff  filter  frequency. 

Table  III  in  the  Appendix  was  constructed  to  facilitate  the 
adjustment  of  R,  C,  and  At  or  the  clock  pulse  frequency.  To  extend  the 
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range  of  possible  clock  frequencies  and  make  it  easier  to  meet  problem 

conditions,  the  built-in  computer  clock  was  passed  through  a  frequency 

divider  with  outputs  equal  to  one,  one-half,  one-fourth  and  one-eighth 
times  the  input. 

If  the  narrow-band  filter  is  in  use,  the  shift  pulse  is  set  to 

make  At  much  smaller,  and  25At  greater,  than  the  inverse  of  the  upper 
frequency, 

F.   INTERRUPT  GENERATOR 

This  was  one  part  of  the  system  treated  with  very  close  attention. 
Since  the  digital  computer  program  may  be  destroyed  by  interrupts  arriving 
at  the  wrong  times,  they  had  to  be  carefully  inhibited  except  when  they 
should  be  expected. 

Figure  3.8  shows  the  interrupt  generator.  As  a  general  rule,  the 
pulse  that  causes  the  interrupt  is  placed  at  the  input  of  a  NAND  gate 
that  has  also  an  input  connected  to  a  line  that  is  kept  off  by  the 
program  except  during  the  sequence  when  it  is  legal  for  the  occurrence 
of  that  particular  interrupt.  In  this  situation  the  output  of  the  gate, 
if  a  pulse  arrives,  actuates  a  delay-flop  that  effectively  generates  the 
interrupt.  Delay-flops  were  used  to  assure  a  short,  well-shaped  pulse. 

Another  interrupt  generator  is  the  circuit  formed  by  the  lines  T013 
(set),  Til 3  [test)  and  the  switch  to  ground.  In  the  subroutine  FASTCOR 
the  digital  computer  tests  line  T113  once  in  eyery   recursion.  If  the 
operator  closes  the  switch  at  the  display  console,  the  line  is  set  low 
and  the  program,  after  reading  this  value,  branches  out  of  the  recursion 
loop,  which  corresponds,  to  a  software  generated  interrupt. 
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G.  AUXILIARY  CIRCUITS 

In  addition  to  the  functional  circuits  described  above,  logic 
circuits  were  developed  to  insure  the  proper  operation  of  the  total 
system. 

Among  these,  the  synchronizer  (Fig.  3.9)  is  responsible  for  keeping 
the  noise  generator  in  the  same  mode  as  the  other  generators.  This 
circuit  also  prevents  the  occurrence  of  sample  interrupts  when  the 
generators  are  not  in  a  COMPUTE  mode.   In  the  Appendix  are  included  the 
logic  equations  represented  by  the  synchronizer.  Other  circuits  and 
connections  are  shown  in  Figs.  3.10  and  3.11. 

The  real -time/time-scale  controls  were  placed  under  digital  command 
through  the  lines  T030  and  T031 ,  respectively. 

A  four-bit  counter  was  patched  to  work  as  a  divider  for  the  built-in 
clock  to  facilitate  the  rate  selection  for  the  feedback-shift  register. 

As  a  sample  clock,  a  frequency  synthesizer  with  square  wave  output 
was  used.  An  analog  comparator  fed  by  this  square  wave  generates  a 
pulse  starting  at  the  zero  crossing  to  trigger  the  sample  interrupt,  or 
to  function  as  a  timer  for  the  recorder  copy  function,  RPLOT, 

The  eight-channel  recorder  was  connected  to  the  outputs  of  the 
signal  combiner,  two  of  the  sine  generators  and  the  pattern  generator. 
Channels  6,  7,  and  8  receive  the  copied  curve  through  the  DAC  lines. 
Channel  5  displays  a  pulse  whenever  the  third  curve  being  plotted  reaches 
the  threshold  setting;  this  pulse  is  generated  by  the  probability  counter, 

The  last  circuit  built  was  a  multiplier-integrator  for  secondary 
computation  of  correlations  by  the  subroutine  ANACOR. 
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IV.  GRAPHICS  TERMINAL 

The  graphics  computer  provides  the  system  one  of  its  most  important 
features,  that  of  human-machine  interaction. 

The  real-time  operation  enhanced  by  the  natural  form  of  graphical 
presentation  made  it  easy  for  the  operator  to  assimilate  information, 
make  decisions,  and  select  the  next  sequence  using  the  system  control 
capability  provided  by  the  graphics  console. 

Since  the  AGT-10  is  a  general  purpose  computer  it  may  be  used  to 
alter  either  the  data  or  the  form  of  the  presentation  to  optimize  the 
man-machine  interface. 

The  system  was  controlled  from  the  graphics  terminal  by  typing  simple 
messages  on  the  screen  using  the  console  keyboard.  A  number  of  control 
possibilities  are  available.  The  operator  could  select  to  put  the 
analog  computer  in  any  mode,  or  test  its  line  values,  select  signals, 
sample,  take  the  Fourier  transform,  perform  auto-  or  cross-correlation 
and  convolution,  or  determine  power  spectrum.  He  could  also  either 
initiate  sequences  to  evaluate  processing  gain,  ROC,  distance,  or  signal- 
to-noise  ratio  or  he  may  manually  control  every   step  of  these  operations. 
Finally,  he  had  a  choice  to  copy  what  was  presented  on  the  screen  either 
on  the  line  printer  or  on  the  paper  recorder. 

In  the  Appendix  is  given  a  more  detailed  description  of  the  operating 
procedures. 
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V.  SONAR  SYSTEM  APPLICATIONS 

As  a  verification  of  the  system  capabilities  several  problems  have 
been  run  to  compare  performances  of  the  three  processors,  check  the 
effects  of  integration  time  and  bandwidth,  and  to  understand  the  signi- 
ficance of  information  and  joint  information  bandpass. 

A.  PRELIMINARY  CONSIDERATIONS1 

To  reproduce  a  signal  from  its  sampled  representation,  in  an  ideal 
case,  it  is  necessary  to  take  samples  at  least  at  the  Nyquist  rate  of 
two  times  the  highest  frequency  contained  in  the  signal.  In  the  case 
of  a  narrow-band-limited  signal,  when  the  bandwidth  is  much  smaller  than 
the  highest  frequency,  the  rate  becomes  twice  the  bandwidth.  The 
recovery  of  the  signal,  however,  requires  filtering  with  steep  filters 
not  found  in  practice. 

For  the  system  described,  the  only  filtering  action  to  which  the 
samples  are  submitted  is  represented  by  the  manner  used  to  generate  the 
vectors  for  display.  This  action  can  be  compared  with  that  of  a  first- 
order  interpolator.  For  this  reason,  whenever  the  shape  of  the  functions 
was  desired,  a  sampling  rate  was  used  of  several  times  the  highest  fre- 
quency, even  in  the  case  of  bandpass  limited  signals  since  no  bandpass 
filter  was  implemented. 

For  the  time-limited  and  bandpass-limited  functions  the  number  of 
independent  samples  (n)  is  twice  the  time-bandwidth  (BT)  product  as  n=2BT. 


Most  of  the  expressions  mentioned  in  this  section  are  from,  or  follow 
directly  from,  those  presented  in  the  theoretical  treatment  of  this  subject 
by  A.  A.  Gerlach,  The  Theory  and  Applications  of  Statistical  Wave-Period 
Processing,  Cook  Electric,  1967. 
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At  least  this  value  is  required  to  exactly  represent  the  function  and 
half  of  that  if  only  the  energy  content  is  the  object  of  the  calcula- 
tions. When  computing  statistical  measures  of  sampled  functions  these 
numbers  must  be  taken  into  consideration;  in  particular,  due  to  the 
relationship  between  energy  and  autocorrelation,  it  is  expected  that  the 
number  of  independent  samples  be  m  =  BT  for  this  case. 

The  practical  limitation  of  the  bandwidth  B  requires  that  when  the 
signal  is  processed  no  noticeable  differences  result  if  the  bandwidth 
is  extended  beyond  that  limit.  This  of  course  depends  on  the  capacity 
of  the  processor  to  "use"  or  "recognize"  the  information  conveyed.   In 
the  case  of  a  matched  filter  or  linear  correlator  this  optimum  bandwidth 
approaches  the  information  bandpass  W  [Ref.  4]  where: 


Ws 


[/   G(f)df]2 


/   G(fTdf 


In  the  general  case  of  detection,  not  only  the  desired  signal  is 
present,  but  also  the  undesired  noise,  and  the  purpose  is  then  to  identify 
signal  in  the  background  of  noise.  The  effect  of  noise  is  to  destroy  part 
of  the  identifiable  information  or  reduce  the  information  bandwidth.  This 
action  is  taken  into  account  when  the  significant  bandwidth  is  taken  as 
the  joint  information  bandpass  or 

OO  CO 

/  G(f)df  /  N(f)df 

w    = ~ 

s,n     °° 

/  G(f)  N(f)  df 

-co 

9 

From  Shannon's  expression  for  the  information  content,  i.e., 
m  =  BT  log2  (1+SNR)  =  [B  log2(l+SNR)  ]T 
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this  interpretation  can  clearly  be  verified.  Therefore,  the  information 
content  of  a  signal  in  the  presence  of  noise  as  seen  by  a  cross- 
correlator  type  detector  may  be  taken  as: 

m  =  WT     where  W  is  the  joint  information  bandpass  for  that 

particular  processor,  signal  and  noise,  and  the  total 
number  of  independent  samples   n  =  2WT. 

It  has  been  determined  that  for  both  the  linear  and  the  amplitude- 
limited  cross-correlator,  W  =  W    ;  and  for  the  clipper-limited 
processor  W  =  1.5  W   .  In  any  case  the  computation  of  W    is  made 
with  the  signals  as  they  are  actually  presented  to  the  correlator,  that 
is,  undistorted,  limited  and  clipped  respectively. 

B,  PERFORMANCE  COMPARISON 

To  compare  the  performance  of  the  three  processors,  an  FM  slide 

signal  was  chosen  with  the  following  characteristics: 

central  frequency  f  =  5000  Hz 

FM-modulation  —  linear  excursion  of  500  Hz  centered  at  f 

sampling  frequency  —  f  '=  1000  Hz 

number  of  samples  —  N  =  256 

pulse  duration  or  integration  time  —  T  =  N/fs  =  256  msec 

Since  W  was  unknown  the  number  of  samples  (N)  could  be  greater,  equal 
or  smaller  than  the  number  of  independent  samples  (n).  However,  consid- 
ering that  the  signal  power  spectrum  is  entirely  contained  in  the  noise 
spectrum  and  assuming  the  noise-spectral -density  constant  over  the  physical 
bandwidth  B,  W  =  B  for  the  linear  case  and  is  greater  for  the  nonlinear 
case.  Thus, 

n  =  2BT  =  2x500x.256  =  256 
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This  was  then  regarded  as  the  linear  processor  having  all  the 
independent  samples  it  was  able  to  process,  while  the  other  two  had 
less  information  than  they  were  capable  of  processing. 

Analysis  of  the  results  showed  that  the  linear  cross-correlator  is 
always  better  than  the  nonlinear  processors  but  it  is  necessary  to 
consider  that  the  last  two  were  processing  less  information  than  they 
could  for  the  preset  bandwidth. 

Also,  it  was  observed  that  for  small  signal-to-noise  ratio  the 
performance  of  the  nonlinear  processors  approached  the  linear  processor. 
This  is  explained  by  inspection  of  Shannon's  expression.  As  SNR  increases 
the  nonlinear  processors  are  unable  to  process  the  amount  of  information 
available  since  the  amplitude  limiting  allows  only  two  identifiable 
states.  This  process  imposes  an  upper  bound  of  2BT  bits  of  information 
as  seen  by  the  processor.  In  other  terms,  since 

SNR  =  f(p,w)SNR,    then 


ASNR  =  ■A.f.(P.?w),AW  SNRT  +  f(p,w)ASNR. 

O      AW  1 


For  the  nonlinear  processors,  when  W  increases  f(p,w)  approaches 
f(p)  =  const  and  ASNR  =  ASNR,  x  const 

while  with  the  linear  correlator  f(p,w)  increases  with  W. 

At  the  lower  end,  if  N  is  such  that  the  full  information  may  be 
processed  by  the  nonlinear  correlators,  the  wider  bandwidth  enhances 
their  performance  and  no  gain  is  achieved  by  the  linear  processor,  since 
the  extra  samples  are  redundant  or  beyond  its  processing  capability. 

C.   INTEGRATION  TIME  X  BANDWIDTH 

In  this  test  the  same  FM  signal  was  used  with  a  reduced  integration 
time. 
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To  decrease  the  pulse  duration,  the  number  of  samples  was  reduced 
to  N  =  128  while  the  sampling  frequency  was  unchanged.  Then  T  became 
128  ms  and  the  information  content  was  cut  to  half  of  the  previous  value. 

Only  the  clipper-limited  processor  was  used.  Comparison  of  these 
results  with  the  others  previously  obtained  revealed  that  the  performance, 
or  processor,  gain  as  well  as  the  detection  sensitivity  was  rediced  as 
expected. 

A  similar  result  could  be  achieved  if  the  pulse  duration  was  kept 
unchanged  and  the  bandwidth  divided  by  two  since  in  this  case  the 
information  content  would  also  be  halved. 

D.   INFORMATION  BANDWIDTH 

In  order  to  explore  further  the  meaning  of  the  information  bandwidth, 
a  CW  burst  code  at  the  frequency  of  5000  Hz  was  generated.  The  number 
of  samples  was  again  256  with  a  sampling  frequency  of  10  kHz  with  a 
resultant  burst  of  25.6  msec.  The  bandwidth  of  this  pulse  was  approxi- 
mately 40  Hz.  The  input  filter  was  kept  adjusted  for  a  bandpass  of  400 
Hz. 

A  second  run  was  also  made  with  a  sampling  frequency  of  100  kHz, 
reducing  T  to  2.56  msec  and  increasing  B  to  400  Hz. 

The  conclusion  derived  from  the  analysis  shows  taht  in  the  first  case 
the  processor  "gain"  was  much  higher  than  in  the  second.  However,  the 
ROC  curves  show  that  this  contributed  nothing  to  the  detection  sensitivity, 
but,  on  the  contrary,  decreased  performance  as  is  shown  by  the  increase  of 
false  alarm  rate.  The  reason  for  this  rests  on  the  fact  that  the  bandwidth 
used  was  much  higher  than  the  information  bandpass,  and  therefore  unneces- 
sary noise  power  was  admitted  over  a  band  where  the  signal  presented  no 
significant  energy. 
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This  confirms  the  concept  that  processor  gain  is  a  significant 
parameter  only  if  the  bandwidth  utilized  corresponds  to  that  compatible 
to  the  processor,  and  to  signal  and  noise  characteristics. 
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VI.  CONCLUSION  AND  RECOMMENDATIONS 

The  final  result  of  the  study  revealed  the  system  as  a  powerful 
instrument  for  signal  analysis.  It  can  be  used  as  an  aid  to  design  work 
or  as  a  complement  for  theoretical  studies  in  signal  processing  and 
detection. 

The  most  important  feature  achieved,  out  of  those  initially  specified, 
was  that  of  the  interaction  between  man  and  computer.  Although  still 
limited  to  some  procedures,  the  operator  can  develop  freely  his  reasoning 
in  a  successive  stream  of  requests.  Signal  analysis  can  then  be  performed 
in  a  conversational  and  profitable  manner. 

In  order  to  have  this  capability,  the  flexibility  of  the  general 
purpose  computers  was  efficiently  explored.  And  no  limits  exist  in  this 
respect.  Many  more  situations  can  be  idealized  and  simulated,  and  other 
applications  devised. 

One  field  that  immediately  deserves  closer  attention  is  that  of  actual 
signal  analysis.  The  system  can  handle  this  task  without  difficulty; 
however,  only  a  continued  experience  can  show  which  improvements  are  to 
be  made.  Depending  on  the  nature  of  the  problem  it  may  be  advisable  to 
make  some  changes  in  order  to  trade  flexibility  and  versatility  for  speed 
and  precision.  In  this  case,  for  instance,  simplified  versions  of  the 
subroutines  FASTCOR,  ANALYSIS  and  PROCESS  could  be  used  without  overlay 
and  with  higher  sampling  rates  to  obtain  very  good  statistical  averages 
in  a  reasonable  time. 

Some  other  applications  the  writer  had  in  mind  to  verify  time  permitting, 
are  commented  upon  next.  They  are  a  suggestion  for  future  studies  with 
the  system. 
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Input  filter:  More  elaborate  filters  either  analog  or  digital  may  be 
simulated  to  represent  a  transducer  accurately. 

Transmission  medium;  The  effect  of  the  transmission  medium  on  the  wave 
propagation  may  be  well  represented  when  the  characteristics  are  known. 
For  instance,  the  frequency  selective  transmission  loss  may  be  repre- 
sented by  a  20  db  per  decade  type  filter  and  the  temperature  and 

salinity  influence  on  the  sound  velocity  could  be  implemented  based 

o 
on  the  equations  giving  c  =  f(T,T  ,S)  [Ref.  9]. 

Multipath  effects  can  be  analyzed  by  using  as  a  signal  a  combination  of 
samples  generated  according  to  a  predetermined  delay  pattern. 

Doppler:  A  very   simple  way  to  simulate  doppler  shift  situations  is  to 
detune  the  pattern  generator  during  the  replica  sampling  periods. 

Statistical  Wave  Period  Processors:  This  processor  simulation  can  be 
obtained  by  a  modification  in  the  subroutine  FASTCOR  to  count  zero 
crossings  of  the  cross-correlation  function  and  then  compute  statistics 
upon  this  data. 

Filter  analysis :  An  example  of  this  possibility  was  worked  with  the 
narrow  band  filter.  Other  filters  either  analog  or  digital  may  be 
analyzed  with  no  difficulty.  If  a  simple  program  is  added,  it  would 
be  possible  to  sketch  desired  filter  waveforms  on  the  screen  and  then 
analyze  them. 

Propeller  signature  recognition:  To  treat  this  very   exciting  problem 
the  harmonic  generator  may  simulate  a  signal  with  a  well-defined  line 
spectrum.  It  is  not  however  an  easy  signal  to  detect  when  it  is 
buried  in  noise,  with  low  signal-to-noise  ratios.  The  recursive  mode 
of  the  correlation  operation  is  appropriate  to  perform  this  computation 
without  occupying  excessive  storage.  This  should  be  a  good  starting 
point  to  develop  a  special  processor  to  handle  this  situation. 
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In  conclusion,  it  is  felt  that  the  signal  processing  system  presented 
here  can  be  an  effective  tool  for  the  exploration  of  processing  methods 
in  both  instruction  and  research. 
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APPENDIX  A 
SYSTEM  OPERATION 

This  Appendix  includes  a  description  of  how  to  control  the  system 
from  the  graphics  terminal. 

Before  loading  the  program,  it  is  required  to  have  accomplished 
the  items  below: 

a)  Assign  a  unit  to  file  3. 

b)  Connect  the  external  clock  (synthesizer)  to  C06  in  the  analog 
board  and  start  the  built-in  clock. 

c)  Connect  the  long  cable  switch  between  the  upper  leftmost  tie 
point  in  the  logic  board  and  ground. 

d)  Load  and  execute  GATED  in  the  AGT. 

After  loading,  the  teletype  prints  a  reminder  to  have  the  analog 
ready,  AGT  GATED  and  clock  running  and  asks  for  ID,  the  terminal 
identification  number.  When  the  message  ID=l(or  2)  is  entered,  the 
operation  starts. 

At  the  terminal  only  the  keyboard  is  necessary  to  type  the  commands. 

In  the  beginning,  the  constants  are  displayed  and  can  be  modified 
by  erasing  the  block  and  typing  the  new  value  with  the  first  digit 
aligned  with  the  first  letter  of  the  heading.  After  carriage  return  the 
SETPOT  mode  is  entered  to  initialize  potentiometers.  Then,  the  system 
waits  for  the  first  command. 

Whenever  an  illegal  operation  or  value  is  typed,  the  message 
"Incorrect  Selection,  Try  Again"  is  displayed  for  a  few  seconds. 

The  commands  are  as  follows: 
A.  CONSTANTS 

M  -  Sample  block  size;  initial  256;  may  be  any  power  of  two  less 
than  or  equal  to  256  and  greater  than  8. 


108 


SAMP  -  1,  2,  or  3  indicates  to  sample  replica,  signal  initial,  or 

signal  continuation;  initialized  with  1. 
THRS  -  threshold  setting;  1000  means  1.0  normalized  threshold;  greater 

values  are  illegal;  initialized  with  0. 
LINE  -  line  number  of  the  sampling  point;  may  be  any  ADK  legal  calling 

number;  initialized  with  1. 
PR0C  -  type  of  processor;  0  is  linear  processor,  4  is  clipper  limited 

and  8  is  amplitude  limited;  initialized  with  0. 
POTS  -  initial  value  of  pattern  generator  injection  potentiometer; 

initialized  with  9999;  the  noise  generator  potentiometer  value 

is  the  complement  to  10000;  in  a  manual  operation  needs  to  be 

set  only  if  the  required  sequence  is  correlation  or  convolution 

recursive. 
STEP  -  value  of  decreasing  step  for  signal  injection  potentiometer  to 

be  used  in  performance  or  ROC  determination;  must  be  less  than 

POTS;  initialized  to  4000. 
CPDT  -  if  5  indicates  that  false  target  rate  is  required;  if  4  , 

initiates  a  distance  determination;  if  3  a  probability  of 

detection;  these  operations,  however  take  place  only  if  a 

correlation  recursive  is  asked  for  next;  initialized  to  one. 
LNTH  -  length  parameter  for  LPL0T;  has  to  be  an  even  number  between 

4  and  16;  initialized  to  16  (smaller  plot  size). 
All  these  are  four-digit  numbers  and  modifiable.  The  next  three  are 
displayed  only. 

SP0W  -  input  signal  power;  initial  1.0 
NP0W  -  input  noise  power;  initial  0 
SNR  (DB)  -  signal  to  noise  ratio  in  db,  initialized  to  0. 
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B.  ANALOG  CONTROL 

SETPOT  -  displays  alternately  the  words  "POT  VAL";  allowable 
potentiometer  identifiers  and  values  may  be  typed  left  justified  under 
the  first  letter;  sets  the  potentiometers. 

SETLINES  -  displays  alternately  "SET  VAL";  allowable  line  identifiers 
(two  digits)  and  1  or  -1  may  be  typed  left  justified;  sets  the  lines. 

TEST  -  displays  alternately  "TST  VAL";  allowable  line  identifiers 
may  be  typed;  after  return  displays  1  or  -1  according  to  the  test; 
requires  another  return. 

POTSET,  COMPUTE,  HOLD  -  puts  the  analog  in  the  mode  specified. 

SCAN  -  displays  current  values  of  all  amplifiers,  potentiometers  and 
several  other  active  components,  read  in  the  current  analog  mode;  waits 
for  the  message  RETN  to  return  or  POTS  to  scan  the  potentiometers  in 
POTSET  mode;  if  the  last  option  was  entered  only  the  potentiometer  values 
are  displayed  and  RETN  must  be  typed  to  return. 

C.  ANALYSIS  MODE 

SAMP  -  asks  for  a  confirmation  about  the  constant  values;  if  NO  is 
typed,  it  brings  in  CONSTANTS;  if  return  is  keyed,  executes  sampling. 
(If  a  recursive  correlation  or  convolution  is  to  be  performed,  the 
replica  must  be  sampled  first  with  the  potentiometer  set  to  9000  and  next 
the  signal,  with  the  signal  input  potentiometer  adjusted  to  the  value 
read  in  POTS). 

FFTR  -  asks  for  function  to  be  transformed  and  displays  function, 
absolute  value  and  phase  of  the  transform;  results  are  not  stored. 

CORRELATION  -  asks  for  type  and  displays  replica,  signal  and 
correlation;  if  recursive,  the  operation  is  terminated  by  the  operator 
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closing  the  manual   switch  for  a  moment;  also  shows  the  threshold  line 
and  information  about  threshold  reached,  probability  of  detection 
(meaningful   only  if  CPDT  were  previously  set  to  3  or  5),  signal-to-noise 
ratio  (the  last  one  computed)   and  distances   (if  CPDT  is  set  to  4). 

CONVOLUTION  -  asks  function  to  be  convoluted  with  the  replica,  and 
type;  displays  replica  flipped,  the  other  function  and  convolution;  if 
recursive,  stopped  by  the  manual  switch. 

POWER  -  asks  for  function  and  type   (.recursive  not  allows);  displays 
autocorrelation  and  power  spectrum  (absolute  value  and  phase);   results 
are  not  stored. 

AUTOCORRELATION  -  asks  for  function  and  displays  function  (twice) 
and  result. 

DISPLAY  -  asks  for  function  and  displays. 

SNRI  -  external  clock  must  be  set  for  time  scale  operation;  result 
may  be  obtained  through  CONSTANTS. 

D.      SIMULATION  MODE 

CWBURST  -  adjusts  the  circuits  for  a  preset  CW  burst  code 

FMSLIDE  -  adjusts  the  circuits  for  a  preset  linear  FM  slide  code. 

FSK  -  sam  for  FSK  code. 

PROPELLER  -  same  for  a  propeller  signature;  no  automatic  operations 
were  implemented  for  this  signal. 

PERFORMANCE  -  computes  the  input-output  characteristics   for  the 
processor  and  signal   selected,  and  displays  the  result;  external   clock 
must  be  set  for  time-scale;  an  extra  return  required  to  end  operation. 

ROC  -  same  as  above  but  with  the  computation  of  Receiver  Operating 
Characteristic 

DIST  -  presets  the  circuits  and  initiates  distance  measurement  operation 
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F.  PLOT 

VPLOT  -  plots  a  VPLOT  on  the  line  printer. 

LPLOT  -  plots  a  LONGPLOT  on  the  line  printer  with  dimension 
inversely  proportional  to  LNTH. 

RPLOT  -  asks  to  ready  the  recorder  and  informs  the  setting  for 
the  probability  counter  potentiometer  if  necessary;  after  return,  informs 
that  a  delay  of  20  seconds  is  in  course  to  allow  starting  the  eight- 
channel  recorder;  during  this  period  a  zero  is  sent  to  channels  6,  7 
and  8  for  calibration;  the  amplitude  scaling  is  to  be  2v  per  division 
and  the  speed  5  or  10  mm.  per  second. 

G.  STOP 

Ends  the  program  and  releases  the  digital  computer. 
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APPENDIX  B 
ANALOG  COMPUTER  INFORMATION 

In  this  Appendix  are  derived  the  logic  equations  to  control  the 
PATTERN  GENERATOR  and  the  FSK  SYNCHRONIZER.  Also  are  included  tables 

for  quick  scaling  the  potentiometers,  (TABLE  I),  for  determination  of 

the  coefficients  of  the  frequencies  generated  by  the  HARMONIC  GENERATOR 

(TABLE  II)  and  for  adjustment  of  the  shift  frequency  of  the  FSR  and  the 
low  pass  RC  filter  (TABLE  III). 

The  following  variables  are  defined: 

RST  -  general  reset  mode 

CMP  -  general  compute  mode 

HLD  -  general  hold  mode 

RRP  -  Ramp  reset  mode 

CRP  -  ramp  compute  mode 

HRP  -  ramp  hold  mode 

R06  -  switch  R06  NC 

RIO  -  switch  RIO  NC 

R26  -  switch  R26  NC 

R25  -  R  input  of  A025 

H25  -  H  input  of  A025 

T07  -  line  T007  (program  call  8) 

T10  -  line  T010  (program  call  9) 

Til  -  line  T011  (program  call  10 

T15  -  line  T016  (program  call  15) 

T17  -  line  T017  (program  call  16) 

C04  -  comparator  C04 

C07  -  comparator  C07 

CWB  -  CW  burst 

FM  -  FM  slide  burst 

FSK-  FSK  burst 

PRO-  propeller  simulation 

DIST  -  distance  burst 
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PTS  -  potset 

SHC  -  clock  line  (FSR-shift) 

RLN  -  reset  line  (FSR) 

CDV  -  frequency  divider 

SIN  -  sample  interrupt  enable 

PATTERN  GENERATOR 


R06  =  T17  =  DIST  +  PRO 


R16  =  T16  =  FSK 


RST  =  T07  •  T10 
CMP  =  TOT  •  TTO 


HLD  =  T07  ■  TIO 

The  above  equations  are  controlled  by  the  program  exclusively  through 
CALL  SETLINES  statements. 

R25  =  HRP  +  CRP 
H25  =  RRP  +  CRP 


HRP  =  HLD  • 

C07  • 

CWB  = 

CWB  •  C07 

•  T07 

■  TIO 

CRP  =  CMP  • 

C07  • 

CWB  = 

CWB  •  C07 

•  T07 

•  TIO 

RRP  =  RST  + 

C07  + 

CWB  = 

CWB  +  C07 

+  T07 

•  TIO 

R25  =  CWB  • 

C07  • 
C07  • 

TIO  • 
TIO  = 

T07  +  CWB 

•  C07 

■  TIO  ■ 

T07 

=  CWB  • 

CWB  •  C07 

•  TIO 

H25  =  CWB  +  C07  +  T07  •  TIO  +  CWB  •  C07  ■  T07  •  TIO 


=  CWB  +  C07  +  T07 


However  C07  is  instantaneously 


C07  =  R25  .  H25  =  CWB  •  C07  •  TTO  <  CCWB  +  C07  +  T07] 


=  CWB  •  C07  •  T07  •  TIO 
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and  then 

H25  =  CWB  +  CWB  •  C07  ■  TOT  •  TTO  +  T07 


=  CWB  +  T07  =  CWB  •  T07 
Making  CW  =  Til  also  controlled  by  the  program,  the  equations: 


R25  =  Til  •  C07  •  T10  ;  H25  =  Til  •  T07  and  RIO  =  C04  were  implemented 

in  the  logic  board  as  shown  in  Fig.  3.3.  To  prevent  possible  occurrence 

of  a  race  problem  and  oscillations  a  short  duration  latch  was  included. 

SYNCHRONIZER 


SHC  =  CMP  +  PTS  +  CDV  =  T07  •  T10  +  PTS  +  CDV 


=  T07  +  T10  +  PTS  +  CDV  =  T07  •  T10  •  PTS  •  CDV 


SIN  =  CMP  •  PTS  =  T07  •  T10  •  PTS 


RLN  =  RST  •  POT  =  T07  •  T10  •  POT  =  T07  •  T10  •  PTS 
The  implementation  is  shown  in  Fig.  3.9. 
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TABLE  I 

POTENTIOMETER  SETTING 

Real  time 

C  = 

.1 

Time  scale 

C  = 

.01 

p   20000  fC 
31 

fsc  =  J       F  =  f\c 

w 


1.59 

9999 

1.19 

7480 

.79 

4960 

.39 

2445 

1.57 

9986 

1.17 

7350 

.77 

4840 

.37 

2325 

1.55 

8874 

1.15 

7230 

.75 

4715 

.35 

2198 

1.53 

8861 

1.13 

7200 

.73 

4585 

.33 

2055 

1.51 

9495 

1.11 

6970 

.71 

4455 

.31 

1948 

1.49 

9352 

1.09 

6840 

.69 

4340 

.29 

1821 

1.47 

9235 

1.07 

6720 

.67 

4205 

.27 

1695 

1.45 

9120 

1.05 

6595 

.65 

4085 

.25 

1570 

1.43 

8995 

1.03 

6465 

.63 

3948 

.23 

1443 

1.41 

8860 

1.01 

6345 

.61 

3835 

.21 

1319 

1.39 

8725 

.99 

6220 

.59 

3703 

.19 

1193 

1.37 

8610 

.97 

6095 

.57 

3580 

.17 

1058 

1.35 

8490 

.95 

5965 

.55 

3478 

.15 

0942 

1.33 

8351 

.93 

5840 

.53 

3325 

.13 

0816 

1.31 

8230 

.91 

5720 

.51 

3205 

.11 

0691 

1.29 

8102 

.89 

5590 

.49 

3078 

.09 

0565 

1.27 

7970 

.87 

5460 

.47 

2935 

.07' 

0439 

1.25 

7851 

.85 

5340 

.45 

2825 

.05 

0314 

1.23 

7730 

.83 

5215 

.43 

2700 

.03 

0188 

1.21 

7600 

.81 

5095 

.41 

2575 

.01 

0063 
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e  =  fC^) 


DC 


TABLE  II 
HARMONIC  CONTENT 


Frequency  Coefficient 


-  K2a2/2  +  3K4a4/8  -  s!<6a6/16 

f  K,a  -  3K3a3/4  +  5K5a5/8 

2f  K2a2/2  -  K4a4/2  +  15K6a6/32 

3f  K3a3/4  -  K5a5/16 

4f  K4a4/8  -  3K6a6/16 

6f  K6a6/32 


e  ■  fCfr  f2) 

Order  Frequency  Amp.   due  to  lowest  K 


1 

Kv 

,K3, 

<K5 

frf2 

^a,  ^b 

2 

Kp  : 

.K4, 

'K6 

2^  ,2f  2 

f1±f2 

1/2  K2a2 
K2a 

3 

K3: 

•K5 

3f] ,3f2 
2f1±f2 

1/4  K3a3 
3/4  K3a2b 

4 

K4: 

•K6 

4f1 ,4f2 

3f1±f2 

2f1±f2 

1/8  K4a4 
1/2  K4a3b 

5 

K5 

5f1 ,5f2 
4f  +f 
3^12^ 

1/16  K5a5 
5/16  K5a4b 

6 

K6 

6f 1 ,6f2 
5f1±f2 
4f1±2f2 
3f-,±3f2 

1/32  K6a6 
3/16  K6a5b 
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TABLE  III 
FSR  FREQUENCY  AND  RC  FILTER 


fcs=V 


Built  in  Clock  Divider  *       fcs  R  C            Series/Parallel 

100  8               1.25  .1  2x10  S 

3x1  P 

100  4               2.0  .1  2x10  S 

100  2               3,33  .1  3x1  P 

100  1               5  .1  2x1  P 

1000  8             12.5  .1  2x1  S 

3x.l  P 

1000  4             20  .1  2x1  S 

1000  2             33.3  .1  3x.l  P 

1000  1             50  .1  2x.l  P 

10000  8           125  3x1  S 

10000  4           200  .1  2x10  S 
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COMPUTER  PROGRAM 
FOR 
HYBRID  SONAR  SIMULATION 

MAIN  PROGRAM 

Purpose:  The  main  program  is  used  for  computer  initialization 
and  serves  as  a  bridge  to  keep  the  status  during  overlay  swapping.  It 
also  creates  two  dummy  pulses  as  Replica  and  Signal  for  operative 
checking  after  initialization. 

Parameters  and  variables:  In  the  explanation  of  the  program  the 

meaning  of  a  variable  is  given  only  in  the  subroutine  where  it  is 

effectively  used.  In  the  case  of  the  MAIN  program  the  variables  below 

are  of  interest: 

M  -  basic  sample  size.  Has  to  be  a  power  of  two.  Initially  has  the 
value  256. 

YR  -  array  for  storage  of  Replica,  within  size  M.  Initially 
dimensioned  to  256. 

YS  -  array  for  storage  of  signal,  with  size  3M.  Initially  dimensioned 
to  768. 

YC  -  array  for  storage  of  correlation,  with  size  2M.  Dimensioned 
to  512. 

VAGO  -  buffer,  with  size  5M  and  dimensioned  to  1280. 

IB  -  buffer,  with  size  2M  plus  one  word.  Dimensioned  to  1025. 

MT  -  text  directory. 

MG  -  graph  directory. 

ID  -  graphics  terminal  number. 
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